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Abstract. 

Context. The fraction of binary stars is an important ingredient to interpret globular cluster dynamical evolution and their stellar 
population. 

Aims. We investigate the properties of main-sequence binaries measured in a uniform photometric sample of 59 Galactic glob- 
ular clusters that were observed by H5T WFC/ACS as a part of the Globular Cluster Treasury project. 

Methods. We measured the fraction of binaries and the distribution of mass-ratio as a function of radial location within the clus- 
ter, from the central core to beyond the half-mass radius. We studied the radial distribution of binary stars, and the distribution 
of stellar mass ratios. We investigated monovariate relations between the fraction of binaries and the main parameters of their 
host clusters. 

Results. We found that in nearly all the clusters, the total fraction of binaries is significantly smaller than the fraction of binaries 
in the field, with a few exceptions only. Binary stars are significantly more centrally concentrated than single MS stars in most of 
the clusters studied in this paper. The distribution of the mass ratio is generally flat (for mass-ratio parameter q>0.5). We found 
a significant anti-correlation between the binary fraction in a cluster and its absolute luminosity (mass). Some, less significant 
correlation with the collisional parameter, the central stellar density, and the central velocity dispersion are present. There is no 
statistically significant relation between the binary fraction and other cluster parameters. We confirm the correlation between 
the binary fraction and the fraction of blue stragglers in the cluster. 



1. Introduction 

The knowledge of the binary frequency in Globular Clusters 
(GCs) is of fundamental importance in many astrophysical 
studies. Binaries play an important role in the cluster dynami- 
cal evolution, as they represent an important source of heating. 
They are also important for the interpretation of the stellar pop- 
ulations in GCs. A correct determination of the stellar mass and 
luminosity functions in GCs requires accurate measure of the 
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fraction of binaries. Stellar evolution in a binary system can be 
different from isolated stars in the field. Exotic stellar objects, 
like Blue Stragglers (BSSs), cataclysmic variables, millisecond 
pulsars and low mass X-ray binaries represent late evolutionary 
stages of close binary systems. The determination of the frac- 
tion of binaries plays a fundamental role towards the under- 
standing of the origin and evolution of these peculiar objects. 

There are three main techniques used in literature to mea- 
sure the fraction of binaries in GCs (Hut et al. 1992). The first 
one identifies binaries by measuring their radial velocity vari- 
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ability (e. g. Latham 1996). This method relies on the detection 
of each individual binary system but, due to actual limits in 
sensitivity of spectroscopy, these studies are possible only for 
the brightest GC stars. Moreover, this technique is sensitive to 
binaries with short orbital periods, and the estimated fraction 
of binaries depends on the eccentricity distribution. The sec- 
ond approach is based on the search for photometric variables 
(e. g. Mateo 1996). As in the previous case, it is possible to in- 
fer specific properties of each binary system (like the measure 
of orbital period, mass ratio, orbital inclination). Unfortunately, 
this method is biased towards binaries with short periods and 
large orbital inclination. The estimated fraction of binaries de- 
pends on the assumed distribution of orbital periods, eccentric- 
ity and mass ratio. Both of these techniques have a low discov- 
ery efficiency and are very expensive in terms of telescope time 
because of the necessity to repeat measures in different epochs. 

A third method, based on the analysis of the number of stars 
located on the red side of the MS fiducial line, may represent a 
more efficient approach to measure the fraction of binaries in a 
cluster for several reasons: 

- The availability of a large number (thousands) of stars 
makes it a statistically robust method; 

- It is efficient in terms of observational time: two filters are 
enough for detecting binaries, and repeated measurements 
are not needed. 

- It is sensitive to binaries with any orbital period and incli- 
nation. 

This latter approach has been used by many authors (e. g. 
Aparicio et al. 1990, 1991, Romani & Weinberg 1991, Bolte 
1992, Rubenstein & Baylin 1997, Bellazzini et al. 2002, Clark, 
Sandquist & Bolte 2004, Richer et al. 2004, Zhao & Baylin 
2005, Sollima et al. 2007, 2009, Bedin et al. 2008, Milone et 
al. 2009, 2010, 2011) to study the populations of binaries in 
individual stellar clusters. The relatively small number of clus- 
ters that have been analyzed is a consequence of the intrinsic 
difficulties of the method: 

- High photometric quality is required and high resolution is 
necessary to minimize the fraction of blends in the central 
regions of GCs; 

- Differential reddening (often present) spreads the MS and 
makes it more difficult to distinguish the binary sequence 
from the single-star MS population; 

- An accurate analysis of photometric errors as well as a cor- 
rect estimate of field contamination are necessary to dis- 
tinguish real binaries from bad photometry stars and field 
objects. 

The first study of binaries in a large sample of GCs comes 
from Sollima et al. (2007), who investigated the global prop- 
erties of binaries in 13 low-density GCs. These authors found 
that the total fraction of binaries ranges from 0. 1 to 0.5 in the 
core depending on the cluster, thus confirming the deficiency 
of binaries in GCs compared to the field where more than half 
of stars are in binary systems (Mayor et al. 1992, Dunquennoy 
& Mayor 1991, Fischer & Marcy 1992, Halbwachs et al. 2003, 
Rastegaev et al. 2010, Raghavan et al. 2010). At variance with 



the high fraction of binaries in field sdB stars (Masted et al. 
2001, Napiwotzki et al. 2004), a lack of close binaries among 
GC hot horizontal branch stars (the cluster counterpart of field 
sdBs) has been confirmed by Moni Bidin et al. (2006, 2009). 

Sollima et al. (2010) extended the study of binaries to five 
high-latitude open clusters with ages between ~0.3-4.3 Gyr and 
found that the fraction of binaries are generally larger than in 
GCs and range between ~0.3 and 0.7 in the core. Very high 
binary fractions have been observed also in some young star 
clusters and for pre-main sequence T-Tauri stars, where the to- 
tal binary fraction might be as high as 0.9 (e. g. Prosser et al. 
1994, Petr et al. 1998, McCaughrean 2001, Duchene 1999). 

These findings suggest that the star formation condition, 
as well as the environment, could play a fundamental role on 
the evolution of binary systems. The binary populations in star 
clusters has been investigated in detail, mainly through Monte- 
Carlo and Fokker-Plank simulations (e. g. Giersz & Spurzem 
2000, Fregeau et al. 2003, Ivanova et al. 2005), N-body (e. 
g. Shara & Hurley 2002, Trenti, Heggie, & Hut 2007, Hurley, 
Aarseth, & Shara (2007), Fregeau et al. 2009, Marks, Kroupa, 
& Oh 201 1) and fully analytical computations (Sollima 2008). 

While the evolution of binaries stimulated by interactions 
with cluster stars could play the major role, there are many 
processes that also influence the binary population in stellar 
systems. For instance binary systems can form by tidal-capture 
(e. g. Hut et al. 1992, Kroupa 1995a). Destruction of binaries 
may occur via coalescence of components through encoun- 
ters or tidal dissipation between the components (Hills 1984, 
Kroupa 1995b, Hurley & Shara 2003). Stellar evolutionary 
processes can significantly effect the property of binaries and 
binary-binary interaction can led to collisions and mergers (e. 
g. Fregeau et al. 2004). The comparison of simulation results 
with observed binary fraction is hence a powerful tool to shed 
light on both the cluster and the binaries evolution. 

In this paper, we report the observational results of our 
search for photometric binaries among GCs present in the 
HST Globular Cluster Treasury catalog (Sarajedini et al. 2007, 
Anderson et al. 2008), which is based on HST ACS/WFC data 
We exploited both the homogeneity of this dataset, and the high 
photometric accuracy of the measures to derive the fraction of 
binaries in the densest regions of 59 GCs. We deserve to future 
works any attempt to interpret the empirical findings presented 
in this paper. 

2. Observations and data reduction 

Most of the data used in this paper come from the HST 
ACS/WFC images taken for GO 10775 (PI Sarajedini), an 
HST Treasury project, where a total of 66 GCs were observed 
through the F606W and F814W filters. For 65 of them, the 
database consists in four or five F606W and F814W deep ex- 
posures plus a short exposure in each band. The pipeline used 
for the data reduction allowed us to obtain precise photometry 
from nearly the tip of the red giant branch (RGB) to several 
magnitudes below the main sequence turn-off (MSTO), typi- 
cally reaching ~0.2 m Q . 

The GO 10775 data set as well as the methods used for 
its photometric reduction have been presented and described 
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in papers II and IV of this series (Sarajedini et al. 2007 and 
Anderson et al. 2008). Q 

The uniform and deep photometry offers a database with 
unprecedented quality that made possible a large number of 
studies (see e. g. Sarajedini et al. 2010 and references therein). 

In this paper we study the main sequence binary popula- 
tion in a subset of 59 GCs. We excluded three clusters (Lynga 
7, NGC 6304, and NGC 6717) that are strongly contaminated 
by field stars and for which there exist no archive HST data 
which could allow us to obtain reliable proper motions and sep- 
arate them from cluster members. We also excluded Palomar 
2 because of its high differential reddening, and NGC 5139 
(oj Centauri), and NGC 6715 because of the multiple main 
sequences (Siegel et al. 2007, Bellini et al. 2010 and refer- 
ences therein). The triple MS of NGC 2808 made the binary- 
population extremely complicated and we presented it in a sep- 
arate paper (Milone et al. 201 la). 

In addition, we also used archive HST WFPC2, WFC3 and 
ACS/WFC images from other programs to obtain proper mo- 
tions, when images overlapping the GO 10775 images were 
available. Table 1 summarizes the archive data used in the 
present paper. 

The recipes of Anderson et al. (2008) have been used to 
reduce the archive ACS/WFC data. The WFPC2 data are an- 
alyzed by using the programs and the techniques described in 
Anderson & King (1999, 2000, 2003). We measured star posi- 
tions and fluxes on the WFC3 images with a software mostly 
based on img2xym_WFI (Anderson et al. 2006). Details on this 
program will be given in a stand-alone paper. Star positions and 
fluxes have been corrected for geometric distortion and pixel- 
area using the solutions provided by Bellini & Bedin (2009). 

2.1. Selection of the star sample 

Binaries that are able to survive in the dense environment of a 
GC are so close that even the Hubble Space Telescope HST is 
not able to resolve the single components. For this reason, light 
coming from each star will combine, and the binary system 
will appear as a single point-like source. In this paper we will 
take advantage from this fact to search for binaries by carefully 
studying the region in the CMD where their combined light 
puts them. 

If we consider the two components of a binary system and 
indicate with mi, m.2, F\, and F2 their magnitudes and fluxes, 
the binary will appear as a single object with a magnitude: 

m bin = mi - 2.5 log(l + 

In the case of a binary formed by two MS stars (MS-MS bi- 
nary) the fluxes are related to the two stellar masses (Ati, Mi), 
and its luminosity depends on the mass ratio q = M2/M1 (in 
the following we will assume M2 < Mi, q < 1). The equal- 
mass binaries form a sequence that is almost parallel to the MS, 

1 Due to a partial guiding failure, we only obtained part of the NGC 
5987 data. In this case the dataset consists in three long exposures in 
F814W and five in F606W, while only the F606W short exposure was 
successfully obtained. For this cluster we obtained useful magnitudes 
and colors for stars fainter than the sub giant branch and with masses 
larger then ~0.2 m Q . 



~0.75 magnitudes brighter. When the masses of the two com- 
ponents are different, the binary will appear redder and brighter 
than the primary and populate a CMD region on the red side of 
the MS ridge line (MSRL) but below the equal-mass binary 
line. 
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Fig. 1. Model MS-MS binary sequences with different mass ra- 
tios for NGC 2298. The dashed-dotted line is the MSRL while, 
continuous black lines indicate the sequences of constant q and 
blue lines mark sequences of constant M\. 

In Fig. Q] we used our empirical MSRL and the mass- 
luminosity relations of Dotter et al. (2007) to generate se- 
quences of MS-MS binary systems with different mass ratios. 

An obvious consequence of this analysis is that our capa- 
bility in detecting binaries mainly depends on the photomet- 
ric quality of the data. Distinguishing the binary populations 
in clusters requires high-resolution images and high-precision 
photometry. Not all stars in clusters can be measured equally 
well. Crowding, saturation, and image artifacts such as diffrac- 
tion spikes, bleeding columns, hot pixels, and cosmic rays can 
prevent certain stars from being measured well. The first chal- 
lenge to this project will be to identify which stars can be mea- 
sured well and which are hopeless. 

In addition to the basic stellar positions and photometry, the 
software described in Anderson et al. (2008) calculates several 
useful parameters that will help us reach this goal. The follow- 
ing parameters are provided for every star: 

- The rms of the positions measured in different exposures 
and transformed into a common reference frame ( rmsx and 
rms\) 

- The average residuals of the PSF fit for each star (^F606w 
and <?F8i4w ) 

- The total amount of flux in the 0.5 arcsec aperture from 
neighboring stars relative to the star's own flux (op606w and 

0F814w)- 
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DATE 


MyC YPTIMF 
IN MAT 1 11V1E 


FILT 


UN o 1 RU 1V1C1N 1 




PI 


ARP 2 


May 11 1997 


Ix260s+5x300s 


F814W 


WFPC2 


6701 


Ibata, R. 




May 11 1997 


5x300s+lx350s 


F606W 


WFPC2 


6701 


Ibata, R. 


NGC 104 


Sep 30 2002 - Oct 1 1 2002 


Ixl0s+6xl00s+3xll5s 


F435W 


ACS/WFC 


9281 


Grindlay, G. 




Jul 07 2002 


6x60s+lxl50s 


F475W 


ACS/WFC 


9443 


King, L R. 




Jul 07 2002 


20x60s 


F475W 


ACS/WFC 


9028 


Meurer, J. 


NGC 362 


Dec 04 2003 


4x340s 


F435W 


ACS/WFC 


10005 


Lewin, W. 




Dec 04 2003 


2xll0s+2xl20s 


F625W 


ACS/WFC 


10005 


Lewin, W. 




Sep 30 2005 


3x70s+20x340s 


F435W 


ACS/WFC 


10615 


Anderson, S. 


NGC 5286 


Jul 07 1997 


3xl40s+lxl00s 


F555W 


WFPC2 


6779 


Gebhardt. K. 




Jul 07 1997 


3xl40s+l 


F814W 


WFPC2 


6779 


Gebhardt, K. 


NGC 5927 


May 08 1994 


6x50s+8x600s 


F555W 


WFPC2 


5366 


Zinn, R. 




May 08 1994 


6x70s+8x800s 


F814W 


WFPC2 


5366 


Zinn, R. 




A f\£. A A 'I 

Aug 06 2002 


30s+500s 


F606W 


ACS/WFC 


9453 


Brown, T. 




Aug 06 2002 


15s+340s 


F814W 


ACS/WFC 


9453 


Brown, T. 




a . . „ ^ o ^ a 1 A 

Aug 2o 2010 


jUs+2X4jjs 


F814W 


T T\ I r C AVI TTV""} 

U VlS/WrL.3 


1 1664 


Brown, T. 




Aug 28 2010 


50s+2x665s 


F555W 


UVIS/WFC3 


11664 


Brown, T. 


JNGC 6121 


J un 1 v 200 J 


15x360s 


r / nw 


ALa/WrL 


9578 


Rhodes, J. 


JNOL. 621o 


r. . ... t A "*>aa A 

Jun 14 2004 


4 X 340s 


F435W 


ACS/ WrC 


1 0005 


Lewin, W. 




Tun i a onn/i 

JUn 14 ZUU4 


ZX4US+ZXOUS 


F625W 


Ape rtyep 


1 0005 


Lewin, W. 


FNOL OJJZ 


in i fine 

ividi zy iyyj 


/XI OUS 


rjjjW 


WrrLZ 


JJOO 


Zinn, R. 




Mar 2y iyyj 


6 x 260s 


F814W 


WrrL.2 


5366 


Zinn, R. 


FNOL OjoB 


Tun 7(1 li ll CVX ~>l\ 1 ( 1 

jun ju - jui uj zuiu 


OXfloUS 


rjyuw 


I r V T ^ IXkTTlC "i 
U V Ij/WrLJ 


1 1739 


r lOttO, \j. 


Mpr A "207 


Ann Al OAA/1 Tun "> ^ IflA^ 

AUg Ul ZUU4 - Jun ZO ZUUJ 


t; v i r, , c v t,i a,, 
J X 1 3 S +J X J4US 


F435W 


a pc rnjnc 

ALi/WrL 


10257 


Anderson, J. 


NGC 6441 


Aug 04-08 20 1 


6x880s 


F390W 


UVIS/WFC3 


1 1739 


Piotto, G. 


iVTrr* a i OA 
INGC t)4yo 


Apr ui lyyy 


inn,. Kill 1A. 

ZX 1 1UUS+4X 1 JUOS 


F606W 


WrrLZ 


6572 


Pare see, F. 




Apr 01 1999 


2x1 100s+4x 1300s 


F814W 


WFPC2 


6572 


Pare see, F 


INUL CO J J 


Ann Ail 1 007 

Aug u4 i yy / 


o X 1 4 US 


F'i^'SW 


WFPC2 


6625 


Buonanno, R. 




Aug 04 1997 


9x160s 


F814W 


WFPC2 


6625 


Buonanno, R. 


NGC 6624 


Oct 15 1994 


6x50s+8x600s 


F814W 


WFPC2 


5366 


Zinn, R. 


NGC 6637 


Mar 31 1995 


6x6()s+8x700s 


F814W 


WFPC2 


5366 


Zinn, R. 


NGC 6652 


Set 05 1997 


12x1 60s 


F814W 


WFPC2 


6517 


Chaboyer, B. 


NGC 6656 


Feb 22 1999- Jun 15 1999 


192x260s 


F814W 


WFPC2 


7615 


Sahu, K. 




Feb 22 1999- Jun 15 1999 


72x260s 


F606W 


WFPC2 


7615 


Sahu, K. 


NGC 6681 


May 09 2009 


32x300s 


F450W 


WFPC2 


11988 


Chaboyer, B 


NGC 6838 


May 21 2000 


2x100s 


F439W 


WFPC2 


8118 


Piotto, G. 




May 21 2000 


2x30s 


F555W 


WFPC2 


8118 


Piotto, G. 


TERZAN 7 


Mar 18 1997 


Ix260s+5x300s 


F814W 


WFPC2 


6701 


Ibata, R. 




Mar 18 1997 


5x300s+lx350s 


F606W 


WFPC2 


6701 


Ibata, R. 



Table 1. Description of the HST additional archive data sets used in this paper, other than those from GO-10775. 



True binary stars will be so close to each other as to be in- 
distinguishable from single stars in our images, so the binarity 
has no impact on the above diagnostics. Therefore, it is safe 
to use the above diagnostics to indicate which stars are likely 
measured well and which ones are likely contaminated. As an 
example, in the six panels of Fig|2] we show these parameters 
as a function of the instrumental [j mp606w and mpsww magni- 
tudes, and illustrate the criteria that we have used to select the 
sample of stars with the best photometry for NGC 2298. 

We note a clear trend in the quality fit and the rms param- 
eters as a function of the magnitude. At all magnitudes, there 
are outliers that are likely sources with poorer photometry and 
that need to be removed before any analysis. Because of this, 
we adopted the following procedure to select the best measured 
stars. 

We began by dividing all the stars of each cluster into bins 
of 0.4 magnitude; for each of them, we computed the median 
values of the parameters rmsx,y and ^F606w,F8i4w defined above 
and the 68.27 th element of the percentile distribution (hereafter 
cr). We added to the median of each bin four times cr, and fitted 
these points with a spline to obtain the red lines of Fig. [2] All 



2 As an example, in the closest GC, NGC 6121, 1 AU corresponds 
to -0.5 mas i. e. -0.01 ACS/WFC pixel. 

3 The instrumental magnitude is calculated as -2.5 log(DN), where 
DN is the total number of digital counts above the local sky for the 
considered stars 



stars below the red line have been flagged as 'well-measured' 
according to that diagnostic. 

The parameters op606w and opsi4w do not show a clear trend 
with magnitude. We flagged as 'well-measured' all the stars 
with o F60 6w < 1 and o F 8i4w < 1. 

In Fig.[3]we compare the color-magnitude diagram (CMD) 
of all the measured stars of NGC 2298 (left), the CMD of stars 
that pass all the selection criteria (middle), and the CMD of re- 
jected stars (right). The sample of stars that have been used 
in the analysis that follows includes stars flagged as 'well- 
measured' with respect to all the parameters we used as di- 
agnostics of the photometric quality. 

The photometric catalog by Anderson et al. (2008) also pro- 
vides the rms of the mp606w and mp8i4w magnitude measures 
made in different exposures. However, a star can have a large 
magnitude rms either because of poor photometry or because it 
is a binary system with short period photometric variability. In 
order to avoid the exclusion of this class of binaries, we have 
not used the rms of magnitude measures as diagnostics of the 
photometric quality in the selection of our stellar sample. 

2.2. Artificial-star tests 

Artificial-star (AS) tests played a fundamental role in this anal- 
ysis; they allowed us to determine the completeness level of our 
sample, and to measure the fraction of chance-superposition 
"binaries". The GC Treasury reduction products (see Anderson 
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Fig. 2. Diagnostic parameters that we have used to select 
the sample of NGC 2298 stars with high-quality photometry. 
The parameters are plotted as a function of the instrumental 
'«F606w and mp8i4w magnitudes. Red circles indicate the me- 
dian rmsx,Y, and <7f606w,f814w per intervals of 0.4 magnitude. 
Red lines separate the well measured stars (thin points) from 
those that are more likely to have poorer photometry (thick 
points). See text for details. 
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sures the images with the same procedure used for real stars and 
produces the same output parameters as in Sect. [2] We consid- 
ered an artificial star as recovered when the input and the output 
fluxes differ by less than 0.75 magnitudes and the positions by 
less than 0.5 pixel. We applied to the recovered ASs the same 
criteria of selection described in Sect. [2]for real stars and based 
on the rms in position and on the <7f606w,f814w and of606w,fsi4w 
parameters. In what follows, including the completeness mea- 
sure, we used only the sample of ASs that passed all the criteria 
of selection. 

Since completeness depends on crowding as well as on stel- 
lar luminosity, we measured it applying a procedure that ac- 
counts for both the stellar magnitude and the distance from the 
cluster center. We divided the ACS field into 5 concentric an- 
nuli and, within each of them, we examined AS results in 9 
magnitude bins, in the interval -14 < mpsnw < -5. For each 
of these 9x5 grid points we calculated the completeness as the 
ratio of recovered to added stars within that range of radius and 
magnitude. Finally, we interpolated the grid points and derived 
the completeness value associated with each star. This grid al- 
lowed us to estimate the completeness associated to any star at 
any position within the cluster. Results are shown in Fig.|4]for 
NGC 2298. The stars used to measure the binary fraction have 
all completeness larger than 0.50. 

3. Photometric zero point variations 

In some clusters, the distribution of foreground dust can be 
patchy, which causes a variation of the reddening with position 
in the field, resulting in a non-intrinsic broadening of the stel- 
lar sequences on the CMDs. In addition to these spreads, small 
unmodelable PSF variations, mainly due to focus changes, can 
introduce slight shifts in the photometric zero point as a func- 
tion of the star location in the chip (see Anderson et al. 2008 
for details). The color variation due to inaccuracies in the PSF 
model is usually -0.005 (Anderson et al. 2008, 2009, Milone et 
al. 2010). In some clusters, differential reddening effects may 
be much larger. An appropriate correction for these effects is 
a fundamental step, as it can greatly sharpen the MS, with a 
consequent improved analysis of the MS binary fraction. 



o o.5 3.1 . Differential reddening 



Fig. 3. CMD of all the measured stars (left), of stars that passed 
our criteria of selection (middle), and CMD of rejected stars 
(right) 

et al. 2008) also contain a set of AS tests. The artificial stars 
were inserted with a flat luminosity function in F606W and 
with colors that lie along the MSRL for each cluster. Typically, 
10 5 stars were added for each cluster, with a spatial density that 
was flat within the core, and declined as r _1 outside of the core. 
The stars were added one at a time, and as such they will never 
interfere with each other. 

Each star in the input AS catalog is added to each image 
with the appropriate position and flux. The AS routine mea- 



In order to correct for differential reddening, we started by 
defining a photometric reference frame where the abscissa is 
parallel to the reddening line, as shown in Fig.[5]for NGC 2298. 
To do this, we have first arbitrarily defined a point (O), near the 
MSTO in the CMD of Panel a. Then we have translated the 
CMD such that the origin of the new reference frame corre- 
sponds to O. Finally, we have rotated the CMD counterclock- 
wise by an angle: 



6 = arctan 



Af6Q6W 



-4f606W-Af814W 



as shown in Fig. [5t>. The two quantities Ap606w and Apsuw 
are the absorption coefficients in the F606W and F814W ACS 
bands corresponding to the average reddening for each GC. 
They are derived by assuming, for each GC, the average E(B - 
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Fig. 4. Left: Completeness as a function of the rapsww magnitude in five annuli (the inner and outer radius of each annulus, in 
ACS pixels, are quoted in the inset) for NGC 2298. Right: Completeness contours in the radial distance versus mpsuw magnitude 
plane. The completeness levels coiTesponding to the red and gray continuous lines are quoted in the figure. Dotted lines indicate 
differences of completeness of 0.05 ranging from 0.05 to 0.95. 
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Fig. 5. Panel a: observed CMD of NGC 2298; the arrow indicates the direction of reddening. The continuous lines are the axes 
('abscissa' and 'ordinate') of the reference frame introduced in the procedure for the measurements of reddening variations. The 
position of NGC 2298 stars in this reference frame is shown in panel b where we draw the fiducial line of the MS as a dashed 
red line. Stars between the dotted lines (black points) have been used as reference stars. Panel C shows the rectified 'ordinate' 
vs. A'abscissa' diagram. 

V) listed in the Harris (1996, 2003) catalog and linearly inter- rotating the CMD is that it is much more intuitive to determine 
polating among the reddening and the absorption values given a reddening difference on the horizontal axis rather than along 
in Table 3 of Bedin et al. (2005) for a cool star. The reason for the oblique reddening line. 
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Fig. 6. Visualization of the local approach for the estimate of the differential reddening suffered by the target star indicated with 
the blue asterisk. Upper left: Position of NGC 2298 stars in the ACS/WFC field of view (gray points). Red circles highlight the 
35 reference stars in the lower left panel. Reference stars are indicated with black crosses, among them, the 35 closest neighbours 
(of the target star) are marked with red circles. The lower left panel is a zoom of a 700x700 pixel centered on the target star. 
Upper right: 'ordinate' vs. A 'abscissa' for all the stars in the NGC 2298 field of view. The median A 'abscissa' of the 35 closest 
neighbor is indicated by the continuous red vertical line and corresponds to the differential reddening value suffered by the target 
star. The histogram of the A 'abscissa' distribution of the 35 closest neighbors is shown in the bottom right panel. 



The value of 6 depends weakly on the stellar spectral type, 
but this variation can be ignored for our present purposes. For 
simplicity, in this section, we will indicate as 'abscissa', the 
abscissa of the rotated reference frame, and as 'ordinate', its 
ordinate. 

At this point, we adopt an iterative procedure that involves 
the following four steps: 

1. We generate the red fiducial line shown in Fig.[5j5. In order 
to determine this line, we used only MS stars. We divided 
the sample of these MS reference stars into 'ordinate' bins 
of 0.4 mag. For each bin, we calculated the median 'ab- 
scissa' that has been associated with the median 'ordinate' 
of the stars in the bin. The fiducial has been derived by fit- 
ting these median points with a cubic spline. Here, it is im- 
portant to emphasize that the use of the median allows us to 
minimize the influence of the outliers as contamination by 



binary stars left in the sample, field stars or stars with poor 
photometry. 

2. For each star, we calculated the distance from the fiducial 
line along the reddening direction (A 'abscissa'). In the 
right panel of Fig. [5] we plot 'ordinate' vs. A 'abscissa' 
for NGC 2298. 

3. We selected the sample of stars located in the regions of the 
CMD where the reddening line define a wide angle with the 
fiducial line so that the shift in color and magnitude due to 
differential reddening can be more easily separated from 
the random shift due to photometric errors. These stars are 
used as reference stars to estimate reddening variations as- 
sociated to each star in the CMD and are marked in Fig. [5] 
as heavy black points. 

4. The basic idea of our procedure, which is applied to each 
star (target) individually, is to measure the differential red- 
dening suffered by the target star by using the position in 
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the 'ordinate' vs. 'abscissa' diagram of a local sample of 
reference stars located in a small spatial region around the 
target with respect to the fiducial sequence. 
We must adopt an appropriate size for the comparison re- 
gion in order to obtain the best possible reddening correc- 
tion. The optimal size is a compromise between two com- 
peting needs. On one hand, we want to use the smallest 
possible spatial cells, so that the systematic offset between 
the 'abscissa' and the fiducial ridgeline will be measured 
as accurately as possible for each star's particular location. 
On the other hand, we want to use as many stars as possi- 
ble, in order to reduce the error in the determination of the 
correction factor. 

As a compromise, for each star, we typically selected the 
nearest 30-100 reference stars and calculate the median 
A 'abscissa' that is assumed as the reddening correction for 
that star. In this way, our differential reddening correction 
will be done at higher spatial frequencies in the more popu- 
lated parts of the observed field. In calculating the differen- 
tial reddening suffered by a reference star, we excluded this 
star in the computation of the median A 'abscissa'. As an 
example, in Fig. [6] we illustrate this procedure for a star in 
the NGC 2298 catalog. The position of all the stars within 
the ACS/WFC field of view is shown in the upper-left panel 
where reference stars are represented by black crosses, and 
the remaining stars are indicated with gray points. Our tar- 
get is plotted as a blue asterisk. The 35 closest neighbor- 
ing reference stars are marked with red circles. The lower- 
left panel is a zoom showing the location of the selected 
stars in a 700x700 pixel box centered on the target. The 
positions of the 35 closest neighboring reference stars in 
the 'ordinate' vs. A 'abscissa' plane are shown in the up- 
per right panel, and their histogram distribution is plotted 
in the bottom-right one. Clearly, neighboring stars define a 
narrow sequence with A 'abscissa' — 0.15. Their median A 
'abscissa', which is indicated by the continuous red line, is 
assumed to be the best estimate of the differential reddening 
suffered by the target star. 

After the median A'abscissa' have been subtracted to the 'ab- 
scissa' of each star in the rotated CMD, we obtain an improved 
CMD which has been used to derive a more accurate selection 
of the sample of MS reference stars and derive a more precise 
fiducial line. After step 4, we have a newly corrected CMD. We 
re-run the procedure to see if the fiducial sequence needs to be 
changed (slightly) in response to the adjustments made and it- 
erated. Typically, the procedure converges after about four iter- 
ations. Finally, the corrected 'abscissa' and 'ordinate' are con- 
verted to mp606w and mpsuw magnitudes. 

From star-to-star comparison of the original and the cor- 
rected magnitudes we can estimate star to star variations in 
E(B - V) and derive the reddening map in the direction of 
our target GCs. As an example, in Fig. [7] we divide the field 
of view into 8 horizontal slices and 8 vertical slices and plot 
A E(B - V) as a function of the Y (upper panels) and X co- 

4 The exact number adopted for each cluster depends on the total 
number of reference stars with a larger number of stars used for the 
most populous clusters. 



ordinate (right panels). We have also divided the whole field 
of view into 32x32 boxes of 128x128 ACS/WFC pixels and 
calculated the average A E(B - V) within each of them. The 
resulting reddening map is shown in the lower-left panel where 
each box is represented as a gray square. The levels of gray 
are indicative of the amount of differential reddening as shown 
in the upper right plot. The analysis of the intricate reddening 
structures in our GC fields is beyond the purposes of the present 
work and will be presented in a separate paper (King et al., in 
preparation). 

Fig.[8]shows the CMDs of twelve of the GCs studied in this 
paper including NGC 2298. These are the clusters that revealed 
the largest differential reddening AE(m F (,06w - f«F8i4w) > 0.05. 




_ 'abscissa' ^'abscissa' 

FG06W Jn Fai4w 



Fig. 7. Bottom-left: Map of differential reddening in the NGC 
2298 field of view. The gray levels correspond to the mag- 
nitude of the variation in local reddening as indicated in the 
upper-right panel. We divided the field of view into 8 hori- 
zontal slices and 8 vertical slices. Upper-left and lower-right 
panels plot A E(B - V) as a function of the Y and X coordinate. 



3.2. PSF Variations 

Some GCs have a reddening that is close to zero and therefore 
we expect negligible variations of reddening within their field 
of view. In these cases, we need to apply only a correction for 
the photometric zero point spatial variation due to small, un- 
modelable PSF variations. Usually, these PSF variations affect 
each filter in a different way, so their most evident manifesta- 
tion is a slight shift in the color of the cluster sequence as a 
function of the location in the field (Anderson et al. 2008). For 
this reason, when the average reddening of the cluster (from 
Harris 1996) is lower than 0.10 mag, we did not follow the 
recipes for the correction of differential reddening described in 
the previous section, but corrected our photometry for the ef- 
fects of the variations of the photometric zero point along the 
chip. We used a procedure that slightly differs from the one 
of Sect. 13.11 The only difference from what done in GCs with 
high reddening is that we did not rotated the CMD and so we 
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Fig. 8. CMDs of twelve GCs before (left) and after (right) the correction for differential reddening. For each cluster we give the 
average reddening from the Harris (1996) catalog. 



did not apply the correction along the reddening line, but along 
the color direction. 

The results of this procedure are illustrated in Fig. [9]where 
we compare the original and the correct CMD of NGC 288. 
The improvement in the quality of our CMD is exemplified by 
the comparison in right panels figures that show a zoom of the 
SGB and the upper-MS. 

Other examples of the improvement in the photometry 
coming from this procedure are shown in Fig.[10]where we plot 
the nine GCs studied in this paper for which we measured the 
largest color variations. The average color variations are typi- 
cally around 0.005 mag for each cluster with E(B - V) <0.10 
studied in this paper and never exceed 0.035 mag. 



4. The measure of the fraction of binaries with 
high mass ratio 

Binaries with large mass ratios have a large offset in luminosity 
from the MSRL and are relatively easy to detect. On the con- 
trary, a small mass ratio doesn't pull them very far off of the 
MSRL, making them hard to distinguish from single MS stars. 



Finally, the low signal to noise photometry of faint stars limits 
the range where binaries can be detected and studied. 

In practice, the limited photometric precision and accuracy 
makes impossible the direct measure of the overall population 
of binaries without assuming a specific distribution of mass ra- 
tios /(q). For this reason, in this paper, we followed two dif- 
ferent approaches to study the binary population in our target 
GCs. 

1) We isolated different samples of high mass-ratio binaries (i. 
e. the binary systems with q > 0.5, 0.6 and 0.7). For them, we 
obtained a direct measure of their fraction with respect to the 
total number of MS stars, and studied the properties of each 
group (Sect. HJl. 

2) We determined the total fraction of binaries by assuming a 
given/(q) (SectO. 

In each cluster, we estimated the fractions of high q bi- 
nary stars in the F814W magnitude interval ranging from 
0.75 (m F814 w,bright) to 3.75 (m F 8i4w,faint) magnitudes below the 
MSTO. In this work we used the MSTO magnitudes from 



In the cases of NGC 6388 and NGC 6441 we used a smaller mag- 
nitude interval between 0.75 and 2.25 magnitudes below the MSTO. 



10 



Milone et al.: The Main-Sequence Binary Population in 59 Globular Clusters 




0.5 1 0.6 0.7 0.8 0.9 



F606W "'Fai4W 

Fig. 9. CMDs of NGC 288 before (top) and after (bottom) correction for photometric zero points variations. Right panels show a 
zoom of CMD around the SGB and upper-MS region. 




Fig. 10. Comparison of the CMD of nine GCs studied in this paper before (left) and after (right) the correction for photometric 
zero points variations. 



This exception is due to the fact that, as we will see in Sect. 14. l.TI we 
do not have reliable proper motions to estimate the numbers of faint 
field stars in the CMDs of these two GCs. 
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Marin-Franch et al. (2009), who used our same photometric 
data base. The choice of this magnitude interval represents a 
compromise between the necessity of a large set of stars and 
the need to avoid faint stars to be able to measure the binary 
fraction also in clusters with poorer photometry (because of 
crowding). 

To illustrate our setup, Fig. [JTJ shows the various regions 
we studied in the CMD of NGC 2298 in order to measure the 
fraction of binaries with mass ratio q > 0.5 for this cluster. 
The upper half of the figure displays two regions of the CMD: 
Region A (upper left) and region B (upper right). 

Region A includes all the stars that we can consider to be 
cluster members. It includes: all the single MS stars and the 
MS+MS binaries with a primary star that have »ip8i4w,bright < 
'«F8i4w < '«F8i4w,faint- The green continuous line is the MS 
fiducial line, drawn as described in Sect. [3] To include stars that 
have migrated to the blue due to measuring error, we extend re- 
gion A up to the green dashed line, which is displaced to the 
blue from the MSRL by 3 times the the average color error for 
a star at that magnitude. The red dotted line is the locus of MS- 
MS binaries whose components have equal mass; we set the 
limit of region A by drawing the red dot-dashed line, displaced 
to the red from the dotted line by 3 times the rms color error. 
The upper-right panel of Fig. QTJshows Region B, which is cho- 
sen in such a way that it contains all the binaries with q > 0.5. 
It starts at the locus of binaries with mass ratio, q-0.5, marked 
by the continuous red line and ends at the dotted-dashed red 
line, which is the same line defined in the upper-left panel. 

The lower half of Fig. [JTJ shows where observed stars and 
ASs fall within these two regions. The left-lower panel plots the 
observed stars and the middle panel shows ASs. We note that 
a significant number of ASs fall in region B. Only a fraction 
of them can be explained by photometric errors; in many cases 
two stars fell at positions so close together that a pair of stars 
has blended into a single object, which would simulate a binary. 
Obviously, in the real CMD, regions A and B are also populated 
by field stars, as shown in the right panel for NGC 2298. We 
will explain how the field star CMD is built in Sect. 14.11 

To determine the fraction of binaries with q>0.5 we started 
by measuring the number of stars, corrected for completeness, 
in regions A (N^ EAL ) and B (N REAL ). They are calculated as 

^real = Tn° BS 1 /ci, where N A ^ is the number of stars ob- 
served in the region A (B) and C\ is the completeness (coming 
from AS tests). Then, we evaluated the corresponding numbers 
of artificial stars (N A RT and N ART ) and field stars (Af EIELD and 
Afp IELD ). In the following Sects. |4"TI and 14.21 we will describe 
the methods that we used to estimate Af EIELD and A^ EIELD and 

The fraction of binaries with q>0.5 is calculated as|j 



xj>0.5 
■'bin 



REAL 



"field 



N A 

REAL 



"art 

■N A N A 
"field "art 



(1) 



6 The first term on the right-hand side of the equation gives the 
fraction of cluster stars (both binaries and blends) observed in Region 
B, with respect to the number of cluster stars observed in Region A. 
The second right-hand term is the fraction of blends and is calculated 
as the ratio of the number of ASs in Regions B and A. 



Similarly, we have calculated the fraction of binaries with q 
>0.6 and q>0.1 . To do this it is necessary to move redward the 
left-hand side (red solid line) of Region B, according to what is 
shown in Fig.[TJ 

The error associated to each quantity of eq.fTJis the Poisson 
error and the error on the obtained binary fraction is calculated 
by following the standard errors propagation. Therefore it rep- 
resents a lower limit for the uncertainty of the binary fraction. 
We note that the binary fractions strongly differ from one clus- 
ter to another with /5 >0 ' 5 ranging from ~0.01 to ~0.40. 

In order to analyze the radial distribution of binary stars 
in GCs and provide information useful for dynamical models 
of our target clusters, we have calculated both the total binary 
fraction and the fraction of binaries with q > 0.5 at different 
radial distances from the cluster center. More specifically, we 
defined three different regions: 

- a circle with a radius of one core radius (rc sample); 

- an annulus between the core and the half-mass radius 
(rc-HM sample); 

- a region outside the half-mass radius (r HM sample). 

The values of the core radius and the half-mass radius are 
from the Harris (1996) catalog. It should be noted that, even 
if our data are homogeneous, in the sense that they came from 
the same instrument (ACS/WFC/HST) and have been reduced 
adopting the same techniques, their photometric quality vary 
from cluster to cluster, mainly because of the different stel- 
lar densities (which affects the crowding). For this reason, for 
some GCs that have poor photometry in their central regions, 
we have measured the fraction of binaries only outside a min- 
imum radius ( /?min) where it is possible to distinguish bina- 
ries with q > 0.5 from single MS stars. The adopted values of 
/?min are listed in Table 2. The fractions of binaries with q>0.5, 

q>0.6, q>0.7 C/j^ ' 5 , /bin° 6 anc ^ /bin" 7 ) f° r tne c l usters m our 
sample are listed in Cols. 3, 4, 5 of Table 2, respectively. In col- 
umn 6 there is also our best-estimate of the total binary fraction 
(i. e. the fraction of binaries in the whole range < q < 1) that 
will be estimated in Sect. 15.21 We give both the fractions of bi- 
naries calculated over the ACS/WFC field and those in each of 
the three regions defined above. 

Following these considerations, it was possible to include 
in the rc sample only 43 out of the original 59 GCs. In addition, 
the limited ACS field of view reduced the number of GCs with 
t"c-hm and r Q HM samples to 5 1 and 45 clusters, respectively. 

4.1. Field contamination 

The best ways to quantify foreground/background contamina- 
tion of regions A and B consists in identifying field stars on 
the basis of their proper motion, which usually differs from 
the cluster motion. For several clusters of the sample consid- 
ered in this paper there are previous epoch HST images with a 
sufficiently long temporal baseline and precision to allow the 
measurement of proper motions. We used archive material to 
determine the proper motions of 20 GCs that are critically con- 
taminated by field stars: ARP 2, NGC 104, NGC 362, NGC 
5286, NGC 5927, NGC 6121, NGC 6218, NGC 6352, NGC 
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Fig. 11. Gray areas in the upper panels are the regions A, and B of the NGC 2298 CMD adopted to select all the (single and 
binary) cluster stars (left) and the candidate binaries with q>0.5 (right), in a range of 3 mp8i4w magnitudes. In all panels, the 
MSRL is represented as a green continuous line, while the green dashed line is blue shifted from the MSRL by three times the 
color error. The red continuous line is the locus of MS-MS binaries with mass ratio g=0.5, while the red dotted line is the locus 
of MS-MS binaries whose components have equal mass. The red dashed dotted line is displaced to the red from the dotted line 
by 3 times the color error. Lower panels show the observed CMD of NGC 2298 (left), the artificial stars CMD (middle), and the 
CMD of field stars (right). 



6388, NGC 6441, NGC 6397, NGC 6496, NGC 6535, NGC 
6626, NGC 6637, NGC 6652, NGC 6656, NGC 6681, NGC 
6838, and TERZAN 7. The procedure to measure proper mo- 
tions is outlined in Sect. 14. 1 .Tl 

In order to determine field objects contamination in the 
CMDs of the remaining clusters, we run a program developed 
by Girardi et al. (2005), which uses a model to predict star num- 
bers in any Galactic field. Details of this procedure are given in 
Sect. l4T2l 

4.1 .1 . Proper motions 

Proper motions are measured by comparing the positions of 
stars measured at two or more different epochs. For the major- 
ity of the clusters only two epochs were available and we fol- 
lowed a method that has been widely described in many other 



papers (e. g. see Bedin et al. 2008, Anderson & van der Marel 
2010). In the cases of NGC 104, NGC 362, NGC 5927, NGC 
6121, NGC 6397, and NGC 6656 we used a sample of images 
taken at three or even more different epochs and determined 
proper motions with the procedure given by McLaughlin et al. 
(2006). We refer the interested reader to these paper for a de- 
tailed description. 

Results are shown in Fig. [12] which plots proper motions 
for twenty GCs. We plotted only stars in the F8 14W magnitude 
range indicated by the numbers quoted in the insets^ Since we 
measured proper motions relative to a sample of cluster mem- 
bers, the zero point of the motion is the mean motion of the 
cluster. Therefore, the bulk of stars clustered around the origin 

7 Note that the magnitude range used to create Fig.[l2]is larger than 
that used to estimate the field star contamination which enters into 
eq. 1 for the the binary fraction calculation. 
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of the vector-point diagrams (VPD) consists mostly of cluster 
members, while field stars are distributed over a larger range of 
proper motions. 

Proper motions offer a unique opportunity to estimate the 
number of field stars that populate the regions A and B of the 
CMD. In order to identify field objects, we began to isolate 
stars whose proper motions clearly differ from the cluster mean 
motion by using the procedure that is illustrated in Fig.QjJfor 
NGC 6656 (where cluster and field stars are well separated in 
the VPD), and in Fig. [Blfor NGC 6838 (where the separation 
is less evident). 

In the left panel of Figs. Qj] and [14] we show the CMD for 
all the stars for which proper motions measurements are avail- 
able. The second column of the two figures shows the VPD of 
the stars in four different magnitude intervals. The red circle 
is drawn to identify the stars that have member-like motions. 
In the following, we will indicate as Rql and /?out the VPD 
regions within and outside the red circles. We fixed the radius 
of the circles at 3.25 <x, where <x is the average proper-motion 
dispersion in the two dimensions. If we assume that proper mo- 
tions of cluster stars follow a bivariate Gaussian distribution, 
the circle should include 99.5 % of the members in each mag- 
nitude interval. The third panel shows the CMD of stars with 
cluster-like proper motion, while selected field objects are plot- 
ted on the right panel. 

We emphasize here that, as we will see in detail in the fol- 
lowing, proper motions are used to evaluate the numbers of 
field stars that randomly fall within the CMD regions A and B 
(Np^ LD ) and not to isolate a sample of cluster stars. This ap- 
proach will allow us to determine the binary fraction by means 
of eq.Q]in the whole ACS/WFC field of view and not only in the 
spatial regions covered by multi-epochs images where proper 
motions are available. To determine the values of N^ LD we 
have to account for three factors: 

1. To accurately measure N^ LD we need a correct estimate 
of the fraction of field stars that share cluster-like proper 
motions. 

2. Proper motions are not available for the whole ACS/WFC 
field of view because, usually, there is only a partial overlap 
between the images at different epochs. As a consequence 
of this we need an accurate measurement of the area of the 
overlapping region. 

3. Proper motions may not be available for a fraction of stars 
in the ACS/WFC catalogs even if they are in the over- 
lapping region because these stars are not measured in 
the second-epoch images (that in many cases come from 
WFPC2), because they either are too faint or in a too 
crowded region. 

Specifically the number of field stars in the region A has 
been evaluated as 
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where: 



have measured proper motions, in the magnitude interval j 
(see the rightmost column of Figs.[T3land fT4l) . respectively. 

- «pj C g LD j is the fraction of field objects that share proper mo- 
tions similar to the cluster; 

- ^area is the fraction of the ACS/WFC field of view with 
multi-epoch observations; 

- c; is the completeness of the ACS/WFC catalog calculated 
in Sect. |Z21 

- /p M is a factor that accounts for the availability of proper 
motions (as in point 3 above). 

And the same is done to evaluate the number of field stars in 
the Region B. In the following, we describe the procedure used 
to determine n RcL , Farea, and / pM . 

Field stars with cluster-like proper motions 
The VPDs of Fig. [12] show that almost all the clusters have 
some field stars that share the mean cluster motion. The frac- 
tion of these sources with respect to the cluster stars depends on 
several factors, such as the astrometric quality of the data, the 
temporal baseline, the line of sight, and the motion of the clus- 
ter with respect to the field. Their fraction is almost negligible 
in NGC 6656 and other cases, but makes a significant contri- 
bution to the binary fraction in most of the GCs of Fig. [T2] 
We now describe a method to determine the fraction of field 
stars with cluster-like proper motion in order to accurately in- 
fer -Vp IELD and A^p IELD in equationQ] 

We note that, for the purposes of this paper, we do not need 
to isolate these intruders. It is sufficient to estimate their total 
amount, and, more specifically, the amount of field stars with 
cluster-like motions that populate the CMD region associated 
with MS -MS binaries or MS single stars. 

We independently calculated, for the GCs with reliable 
proper motions, the number of field stars with cluster-like 
proper motions for each of the four magnitude intervals of 
Figs. [T3] and [14] In the cases of GCs where cluster and field 
stars are clearly separated in the proper motion diagram (ARP 
2, NGC 104, NGC 362, NGC 5286, NGC 6121, NGC 6218, 
NGC 6388, NGC 6397, NGC 6496, NGC 6535, NGC 6637, 
NGC 6624, NGC 6652, NGC 6656, and Terzan 7) we used the 
method that is illustrated in Fig.[l5]for NGC 6656. All the field 
and cluster stars with reliable proper motions are located within 
the dotted circle of the left panel VPD We considered as proba- 
ble cluster members all the objects that are plotted as thin gray 
dots in the yellow area (region Rql) of the zoomed VPD in the 
right panel, while remaining objects are flagged as field stars 
and are represented as heavier points. 

The distribution of field stars in the VPD is clearly elon- 
gated and the isodensity contours can be approximately de- 
scribed by ellipses. In Fig. [15] we show the two isodensity con- 
tours that are tangent to the region Rcl and define the red region 
(Rj). The number of field stars within Rql is assumed to be: 

. Rn 



'FIELD 



"R t Sr cl /Sr t 



- «FiELD,j and ^ field j Slle t ^ le tota ^ numrjer of field objects and 
the number of field objects within Region A for which we 



where Sr cl and Sr t are the areas of regions Rcl and Rj and 
«r t is the number of stars within Rj. 

In the cases of NGC 5927, NGC 6352, NGC 644 1 , NGC 668 1 , 
and NGC 6838, where the separation of field and cluster stars is 
less evident, we followed a different recipe, which is illustrated 
in in Fig.[l6]for NGC 6838. The upper panels show the CMD 
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Fig. 12. Vector-point diagram of proper motions, in equatorial coordinates, for twenty GCs. Only stars in the mp8i4w interval 
indicated in each panel are shown. 
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(left) and the VPD {right) for stars in the third interval of mag- 
nitudes (j=3) of Fig. [14] We selected, on the CMD, a sample of 
stars that, on the basis of their color and magnitude, are proba- 
ble background/foreground objects. These stars are marked as 
heavy black points in the lower CMD of Figure[16] while in the 
right-lower panel we show their position in the VPD. 



If we assume that the fraction of selected objects within 
Rcl with respect to the total number of selected field ob- 
jects (/j^eld) * s re P resen t a ti ve of the overall fraction of 
field stars that share cluster proper motions we can impose: 



"held = "field/held- The contribution of n^ LD to the 



measure of the binary fraction is, 
than 0.01. 



for all the clusters smaller 
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Fig. 13. Leftmost column: mpsww vs. mp606w - '«F8i4w CMD for all the stars of NGC 6656 with available measures of proper 
motions. Second Column: Proper motion diagrams of the stars in the left panels in masyr~ l , in intervals of 1.4 mag. Third column: 
The proper motion selected CMD of cluster members. Rightmost column: The CMD of field stars. 
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Fig. 14. As in Fig. [Bj but for NGC 6838. 



In order to investigate the reliability of this approach, we 
applied it also to the 15 GCs for which proper motions allow 
us to almost completely separate cluster stars from field ones. 
In all cases, we found full consistency between the two ap- 
proaches, with the the fraction of binaries with g>0.5 listed in 
Table 2 differing by less than 0.01. 

The spatial coverage of multi-epoch images 

For most clusters, there is only a partial overlap among the dif- 
ferent epoch images. In the following we will refer to the region 
that has been observed in at least two epochs as 'i?n'. Fig. [TTl 
shows the example for NGC 6656, where we indicate as light 
gray points all the stars for which we have only photometry, 
and mark with black points the stars with both photometric and 
proper motion measurements. As our field is just a few square 
arcmins, we can assume that the background/foreground popu- 



lation is uniformly distributed within it, and therefore we esti- 
mate the total number of field stars in our field of view as the 
product of the number of field stars in the region Z?n and the 
ratio between the area of the total field of view and the area of 
Rn, In this paper, we will refer to this ratio as: Farea- 

Completeness correction for field stars 

In the procedure that we have applied to determine the clus- 
ter membership using proper motions, we have automatically 
excluded all the stars that might be members but have poor as- 
trometry. An accurate estimate of the fraction of these stars is 
necessary to infer the correct values of A^p IELD and A^ IELD . To 
estimate the fraction of cluster stars lost by applying the proper 
motion selection criteria, we applied the procedure illustrated 
in Fig. [18] for NGC 6656. In panel a we show the mp8i4w vs. 
»iF606w - ?WF8i4w CMD for all the stars in the region 'Rn. Proper 
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Fig. 16. Estimate of the contamination of field stars that 
share cluster proper motions in NGC 6838. Upper panels 
show the CMD (left) and the VPD (right) for stars with 
18.15<m F8 i4 W <19.45. Lower panels display the CMD (left) 
and the VPD (right) for those objects that, on the basis of their 
position on the CMD, are probable field stars. 
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Fig. 17. Map of all the stars NGC 6656 (light gray points) with 
7«F8i4w < 18.25. Black points mark all the stars with available 
proper motions. The contour of region 'Ru' is colored in red. 

motion measurements are available only for a fraction (/pm) of 
these stars. Their CMD is shown in panel b, while the CMD for 
stars with no available proper motions is plotted in panel (c). 

To determine /pm we started by dividing the CMD into bins 
of 0.5 ;«F8i4w magnitudes. In each of them, we counted the total 
number of observed stars (Nobs) and the number of star with a 



reliable estimate of proper motions (Npm)- The fraction of stars 
with a proper motions in that bin is:/pM = Npm/Nobs- 

We then calculated the median otfsmw magnitude of the 
observed stars (tomed) in each bin. We associated to each bin 
the corresponding value of/pM and ;«med- The (J^, M ) for each 
i-star is calculated by interpolation with a spline. In panel (d) 
of Fig. [18] we show the final /pm as a function of mp8i4w- For 
the GCs studied here always we have/pM > 0.4 at the level of 
3.75 F814W magnitudes below the MSTO. 
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Fig. 18. Panel a: CMD for all the stars within the region 'Ri{ ; 
Panel b: CMD for stars with a reliable estimate of proper mo- 
tions; Panel c: CMD for stars within 7?n' for which there are 
not reliable measurements of proper motions; Panel d: Fraction 
of stars with good proper motions with respect to the total num- 
ber of stars within 'i?n' as a function of m F 8i4w- 



4.1 .2. Galactic model 

In order to estimate the number of background/foreground stars 
in the field of view of the GCs studied in this paper, and for 
which we do not have reliable measurements of proper mo- 
tions, we used the theoretical Galactic model described by 
Girardi et al. (2005). This model was used to generate a syn- 
thetic CMD (in the ACS/WFC F606W and F8 14W bands) con- 
taining the expected field stars in the cluster area that we are 
studying. The synthetic CMDs were used to count the number 
of field stars in the CMD regions A, and B (N^ m , N$ m ) defined 
in Fig Q~TJ Obviously, the number of stars in simulated CMDs 
may differ from that of observed field stars. To minimize the 
effect of such uncertainties on the measure of the fraction of 
binaries in GCs, we defined in the CMD a region F on the red 
side of equal-mass binaries fiducial sequence, that is delimited 
on the blue side by the red dashed-dotted line of Fig. [TTJ and 
is likely not populated by cluster stars, as illustrated in Fig. [19] 
for NGC 2298. We determined the numbers of stars within F 
in the observed and in the simulated CMDs (Nq BS and N^ m 
respectively). 

The number of field stars in the CMD regions A is then 
calculated as: 

^F A IELD=^ S A IM ^ BS /^V S F IM (3) 

and a similar equation is used to estimate the number of field 
stars in the region B. 
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As anticipated in Sect. 2, we removed from our list all clus- 
ters for which we had no proper motion (two epochs data) and 
for which Girardi et al. (2005) model was prediction a field 
star contamination larger than 1 %, with the only exception of 
E3 (a 2.4% expected contamination) and NGC 6144 (1.3%). 
Therefore, for clusters for which we have to rely on a Galactic 
model to estimate the foreground//background stars, the con- 
tamination is expected to be minimal. On the other hand, we 
kept into the sample all cluster for which we could use proper 
motion to estimate field stars, independently from the level of 
contamination. In order to investigate whether the estimate of 




Fig. 19. Observed CMD of NGC 2298 (left) and simulated 
CMD of field stars in the NGC 2298 field of view (right). 



field stars from Galactic models is reliable, we applied the syn- 
thetic CMDs method also in the 15 GCs for which we have reli- 
able proper motion measurements. We found that, in the cases 
of GCs with a small field-star contamination, the fraction of 
binaries with q>Q.5 derived following the two approaches is 
identical within the uncertainties, with differences smaller than 
0.01. For some GCs with a significant background/foreground 
population, namely NGC 5927, NGC 6352, NGC 6388, NGC 
6441, NGC 6637, and NGC 6681, the fractions of binaries de- 
rived using a Galactic model differ from those derived using 
proper motions by 0.01 to 0.03 (for NGC 6441). 

4.2. Estimate of the fraction of apparent binaries 

Chance superpositions of two physically unrelated stars that 
happen to lie nearly along the line of sight (apparent binaries) 
and superposition of a faint star and a positive background fluc- 
tuation may reproduce the color and luminosity of a genuine 
binary system, and populate the CMD region occupied by bi- 
naries. In a crowded stellar field, like the core of a GC, a reli- 
able measure of the binary fraction requires good accuracy in 
deriving the number of chance superpositions. 

We can identify and reject a significant fraction of these 
objects by analyzing the stellar profile, and the PSF-fit errors. 
For this reason, in this work, we limited our study to the objects 
that pass the selection criteria described in Sect. 12.11 

In order to account for the blends that have not been re- 
jected, a statistical estimate of their number and distribution 



in the CMD is necessary. In this paper, we used extensive 
artificial-star test experiments to evaluate directly the effects 
of blends. 

Specifically, in this subsection, we illustrate the procedure 
adopted to determine the relative numbers of artificial stars in 
the regions A and B of the CMD of Fig. EE]J(A^ RT and ^art) 
that are used to calculate the last term of eq.Q] 

This analysis requires that the artificial star sample that we 
will compare to observed data reproduce as much as possible 
all the details of real stars. In particular we need the best possi- 
ble match between the luminosities, the radial distribution and 
the photometric errors of observed and simulated stars. 

The data set described in Anderson et al. (2008) includes 
an extensive set of artificial-star tests for each cluster. The same 
quality parameters were determined for the artificial stars as for 
the real stars, so we apply the same selection criteria to them as 
we did to the real stars in Sect. 12.11 

To apply these generic artificial-star tests to the real clus- 
ter distribution, for each real star observed, we took a set of 
the artificial stars within +0.10 magnitude and with radial dis- 
tances within 100 pixels of that of the star. These are the stars 
that were used to estimate the measurement errors (random and 
systematic) of the stars in the cluster. 

The result of this procedure is a catalog of simulated stars 
that reproduces both the radial and the luminosity distributions 
of real stars. Several effects contribute to the observed width 
of the main sequence. In addition to photon noise, we have the 
contribution of spatial variations of the PSF and residual differ- 
ential reddening that are beyond the sensitivity of the method 
that we used to correct them, as well as scattered light, pos- 
sible star-to star metallicity variations, etc. However, for the 
purposes of this work, it is not necessary to distinguish the con- 
tributions of the single sources of the broadening and we can 
include them in the photometric errors (cr). 

Since MS-MS binary systems and apparent binaries both lie 
on the red side of the MS, we can use the MS scatter to the blue 
side of the MS as an estimate of the photometric error. We note 
that the blue portion of the MS may be contaminated by MS- 
white dwarf binaries but their influence on cr is expected to be 
negligible, and further reduced by the applied "kappa-sigma" 
rejection algorithm, as described below. 

In order to estimate cr, we used the following iterative 
procedure, which has been applied to both the observed and 
artificial-star CMD. First of all, we subtracted the color of the 
MSPvL from the color of each star. Then we divided this CMD 
into several intervals of magnitude, each one containing the 
same number of stars, and constructed a histogram of the color 
distribution for each magnitude interval. The size of each inter- 
val is a compromise between maximizing the number of stars 
to reduce the statistical errors and minimizing the magnitude 
intervals to account for the variations of the photometric error 
as a function of the luminosity. For these reasons, the size of the 
adopted interval varies from one cluster to another, depending 
on the number of sampled stars. 

We used least-squares to fit each histogram with a Gaussian 
that had three fitting parameters: its center, its amplitude, and 
its dispersion cr. Then, we rejected all stars for which color is 
far more than 3cr from the fiducial line, because most of these 
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objects must be field stars or binaries. Finally we used the re- 
maining sample for a new Gaussian fit. 

All the stars with negative color in the rectified CMD (i. e., 
those on the blue side of the MS) are used for a new Gaussian 
fit, but, this time, we fixed the center and the amplitude of the 
Gaussian and considered cr as the only free parameter. The best 
fitting cr is adopted as the average photometric error in that 
magnitude interval. The errors corresponding to a given mag- 
nitude in the CMD are obtained by interpolations. 

As expected, the artificial star color distribution is narrower 
than the real star one. We need to properly estimate the differ- 
ence between the artificial-star photometric error and the pho- 
tometric error of real stars, since, as it will be clearer in next 
section, we need an artificial-star CMD with the correct photo- 
metric error in order to estimate the photometric outliers which 
contaminate the binary region. 

The smaller color dispersion of the artificial star CMD 
comes from the fact that the measurement errors of artificial 
stars are smaller than the corresponding error of real stars. This 
difference is due to the fact that, in fitting artificial stars, we use 
exactly the same PSF that was used to originate them, while 
we cannot expect the same perfect match of the PSF with the 
real PSF of real stars. In addition, and for the same reason, 
artificial-star photometry is not affected by zero point photo- 
metric errors, and errors associated with the differential red- 
dening correction. 

The difference between the MS color spread of observed 
and simulated stars might be also due to multiple stellar pop- 
ulations. Indeed, nearly all the GCs studied so far host two or 
more generations of stars with a different light-elements. In few 
GCs, there are also star-to-stars iron variations (see Milone et 
al. 2010b for a recent review). 

Among the clusters studied in this paper, multiple MSs as- 
sociated to helium variation have been identified in 47 Tuc, 
NGC 6752, and NGC 6397 where the »iF606w - '«F8i4w color 
difference between the He-rich and He-poor MS is about 0.01 
mag (Anderson et al. 2009, Milone et al. 2010a, 2011b,c) i. e. 
has the same order of magnitude as the color errors of the best 
measured MS stars. NGC 6656 (M22) is the only cluster of this 
paper where two groups of stars with a different iron content 
have been identified. In this case theoretical isochrones show 
that the measured [Fe/H] difference of ~0.15 dex do not pro- 
duce any appreciable otf(,06w - »Jf814w color bimodality among 
MS stars (Marino et al. 2009, 2011). In general the MSs cor- 
responding to the different stellar populations observed in the 
majority of GCs (and hence formed by stars that could have dif- 
ferent overall CNO abundance, and light elements variations) 
are almost overimposed when observed in the »Jf606w - '«F8i4w 
color (Sbordone et al. 201 1). 

As an example, the difference in color dispersion between 
the real and the artificial star CMDs of NGC 2298 are shown in 
Fig- EHl In order to compare the real and the artificial star color 
distribution it is necessary to appropriately re-scale the latter. 
For this, we considered the measured dispersions as a func- 
tion of the mp8i4w magnitude for both observed and simulated 
MSs, and calculated by least squares the 4 th order polynomials 
(Preal and Parts) that best fit each of them. As an example, 
Fig- IZU (upper panel) shows the measured dispersions and the 



best fitting functions for the case of NGC 2298. In this paper, 
we considered the spread of the MS stars as a reliable indicator 
of the photometric errors to be associated to color measures. 
We believe that it represents a much more accurate estimate for 
the observed MS breadth than the one given by the rms value 
obtained from magnitude measures of the single AS MS stars. 
In fact it also accounts for residuals photometric zero point er- 
rors, errors associated to the reddening correction method and 
possible intrinsic spread due to the presence of multiple stellar 
populations. 

The difference between the obse rved and simula ted MS dis- 
persion is expressed as: Acr VI = JP^^ - ^arts- Assuming 
that any spread of MS stars around the MS fiducial line comes 
only from photometric errors, Acr VI indicates how the artificial- 
star color errors underestimate our real-star photometric error. 
As a final, fundamental step for the following discussion, we 
made the artificial-star CMD similar to the observed one by 
adding to each artificial star additional random noise in color, 
extracted from a Gaussian distribution with dispersion Acr VI . In 
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Fig. 21. Upper panel: MS dispersion as a function of mpsi4w 
magnitude for NGC 2298. Circles and triangles correspond to 
real and artificial stars, respectively. The corresponding best fit- 
ting 4 th order polynomials are represented by dotted and dashed 
lines. Lower panels: comparison of the observed CMD of NGC 
2298 (left) and the CMDs of artificial stars before (middle) and 
after (right) the application of color spread. 



the bottom panels of Fig.|2T|we compare the observed CMD of 
NGC 2298 and the CMDs of artificial stars before and after the 
application of the color spread. The latter CMD has been used 
to calculate the ratio between A^ RT and A^ RT used in Eq.Q] 
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Fig. 20. Left: The observed (top) and simulated (bottom) CMD of NGC 2298 with the fiducial line overplotted; Middle: The 
CMDs rectified by subtraction of the fiducial line; Right: Color distribution of the rectified CMDs. The <x in the inset are those 
of the best-fitting Gaussian. 



5. Results 

In this Section we illustrate and discuss the main results of this 
work. Specifically: 

- In Sect. 15.11 we analyze the mass-ratio distribution of bi- 
naries in each of the 59 GCs studied in this paper in the 
range 0.5 < q <1. Results from individual clusters are used 
to estimate the average mass-ratio distribution of binaries; 

- Attempt to calculate the total fraction of MS-MS binaries 
is proposed in Sect. 15.21 

- Sect. 15.31 gives a summary of the literature measurements 
of the binary fraction in GCs and compares them with ours. 

- In Sect. 15.41 we investigate the distribution of binaries as a 
function of the primary star mass (magnitude); 

- The radial distribution of binaries in each GC is studied in 
Sect. 1331 

- Finally, monovariate relations between the binary fraction 
and the main parent cluster parameters (absolute luminos- 
ity, central velocity dispersion, metallicity, age, central den- 
sity, ellipticity, core and half mass relaxation time, HB mor- 
phology, collisional parameter) are discussed in Sect. 15.61 



5.1. Mass-ratio distribution 

This section, presents the mass-ratio distribution of the binary 
population for our target GCs in the range of 0.5<q<l. To do 
this, we have divided Region B of the CMD into five inter- 
vals of mass ratio (fi/,2,...,j) as shown in Fig.|22]for NGC 2298. 
We chose the size of these regions in such a way that each of 
them covers almost the same area in the portion of the CMD 
populated by binary systems with q > 0.5. The sub-region Z?j 
includes also the gray area on the right side of equal-mass bi- 
naries fiducial that is populated by binary systems with q ~ 1 
but large photometric errors. 

The fraction of binaries in each sub-region B\ is calculated 
over the entire WFC field of view following the procedures de- 
scribed in Sect. [4] Each sub-region includes binary stars within 
a given mass-ratio interval (Aq\) as labeled in Fig. [22] To ac- 
count for the different mass-ratio values of each sub-region, 
and analyze the mass-ratio distribution, we derived the normal- 
ized fraction of binaries: 

Vbin,i = /bin,i/A<7i.0 

8 If we assume that: 
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Results for all clusters are shown with black symbols in 
Figs. |231 and l24l To increase the statistics, we have also divided 
Region B into two large mass ratio intervals with 0.5 < q < 0.7 
and 0.7 < q < 1 and calculated Vbi n in each of them. The results 
we obtained by using these q bins are marked with red open 
circles in Figs. l23l and l24l 




Fig. 22. As an example, we show with gray and black shaded 
areas the five CMD regions (.61,2,. ..,5) used to determine the 
mass-ratio distribution of binary stars in NGC 2298. 

The mass-ratio distribution is almost fiat for most of the 
GCs of our sample but in few cases we cannot exclude possible 
deviations from this general trend. To investigate this statement 
we compared the observations with a flat distribution, calcu- 

- <f>(q) is the continuous function that describes the distribution of 
the number of binaries as a function of the mass ratio. 

- N is the total number of stars (both binaries and single stars) 

- 2V21, iV5J, N?? the number of binaries in each region Bi 2 ,...,s ■ 

Obviously 

X' 5 <t>(<l)dq=f^ ^dq+J^ (f>(q)dq+...+J^ <p(q)dq 
where [ql:q2], [q2:q3], [q5:q6] are the mass-ratio intervals corre- 
sponding to the CMD regions of Fig. [22] We have: 

At this point, the best we can do to gather information on <f>(q) is to 
use the approximation: 

<f>(q)dq = (/>*(q)(q(i + 1) - q(i)) = (/>*(q)&q, 
and calculate: 

v w = ( C v> = c/a®. 

If we normalize c/>*(q) by the total number of stars we find that 
the normalized fraction of binaries differs from <f>* by a factor 1/N: 
<f>*(q)/N = N^/iNAqd = f^JAqt = v bmj . 

Since the total number of stars changes from one cluster to each other, 
we use here Vbm.i as the best approximation of the mass-ratio distribu- 
tion in each q interval. 



lated for each cluster the reduced^ 2 and quoted it in Figs. [23] 
and [24] Montecarlo simulations demonstrate that in the case 
of a flat distribution we expect the 50% of the total number 
of clusters having x 2 <U1 an d the 99% x 2 <3.8. We found 
X 2 values higher than 3.8 in four GCs namely NGC 6366 ( 
* 2 =4.92), NGC 6496 Cr 2 =6.38), TERZAN 7 (^ 2 =4.45) and 
E 3 Cf 2 =13.62). 

To compare the trend of the fraction of binaries as a func- 
tion of q for different GCs we divided Vbi n ,i by two times the 
fraction of binaries with g>O.50. 

Results are in Fig. [25] Black points indicate the measure- 
ments for all the GCs, while red points with error bars are the 
averages in each mass-ratio bin. The gray line is the best fitting 
line. Its slope is indicated in the figure and suggests that the 
mass-ratio distribution is nearly flat for g>0.5. In the Appendix 
we will demonstrate that this result is not affected by any sig- 
nificant systematic error. 

Since we have determined the mass-ratio distribution over 
the entire ACS/WFC field of view, our conclusions should indi- 
cate the general behavior of the binaries in GCs. Unfortunately, 
due to the relatively small numbers of binaries, we could not 
extend this analysis to each sample of rc, the rc-HM, and the 
>"oHM stars. In these regions, due to mass-segregation effects, the 
mass-ratio distribution could differ from that shown in Fig. [25] 

Up to now, there are few observational constraints on the 
overall mass-ratio distribution of the binary population in GCs. 
One of the few measures of/(q) for binary systems, available 
in the literature, comes from Fisher et al. (2005) who estimated 
the the mass-ratio distribution function from spectroscopic ob- 
servations of field binaries within 100 parsecs from the Sun. 
The/(q) derived by Fisher et al. (2005) is shown in the upper 
panel of Fig. [26] Binaries with q<~0.9 have a nearly flat distri- 
bution while there is a large concentration of binaries formed 
by two components of similar mass. Tout (1991) studied the 
binary systems located in the local field and suggests that/(q) 
can be derived by randomly extracting secondary stars from the 
observed initial mass function (IMF). The mass-ratio distribu- 
tion that we obtain by randomly extracting pairs of stars from 
a Kroupa (2002) IMF is displayed in the upper panel of Fig. [26] 
for MS binaries with a primary with 0.47 < M < O.76M 
which is the typical mass interval corresponding to the mag- 
nitude interval we analyzed in the present work. In this case, 
the /(q) shape rapidly decreases from low to high mass-ratio 
values with only the 24% of binaries having q>0.5. 

In order to investigate whether the observations of Fig. [25] 
are consistent with any of the two mass-ratio distributions de- 
scribed above, we calculated the normalized fraction of bina- 
ries we expect in the CMD of a GCs where binary stars fol- 
low the distribution by Fisher et al. (2005) and the distribu- 
tion obtained from random extraction of secondary stars from 



9 since v bm i depends on the fraction of binaries, which changes from 
one cluster to each other, to compare results from different clusters, we 
have to normalize it by the total fraction of binaries. Due to the lack of 
information on binaries with q<0.5, we normalized v bm i by f£* ' s . We 
also multiplied the latter by a factor of two to normalize to one. (Note 
that, by chance, 2/^ ' 5 corresponds to the total fraction of binaries for 
the case of flat mass-ratio distribution). 
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Fig. 23. Mass-ratio distribution for the binary population in the ACS field of 29 GCs. Black filled-circles show the normalized 
fraction of binaries in five mass-ratio intervals, while red open-circles indicate the Vbin,i values obtained by using only two bins 
with 0.5<g<0.7, and 0.7<q<l. Horizontal segments mark the adopted mass-ratio interval, while observational errors are plotted 
as the vertical lines and shadowed areas. 



a Kroupa (2002) IMF (vbi n ,F, vm^r). ^ e a ^ so divided each of 
these quantity by two times the fraction of binaries with q>0.5 
of the corresponding CMD C/j^rf, fwfaji) m close analogy to 
what done for real stars. 

Results are in Fig. [27] where the values of Vbin,F/(2 /^p 5 ) 

and Vbin,R/(2 f^ R ) are plotted as a function of q. The best- 
fitting least-squares lines are colored gray and their slopes are 
quoted in the inset. Red points are the observed average binary 
frequencies of Fig. [25] The large reduced-^ square values ob- 
tained from the comparison of the theoretical and the observed 
points, and quoted in the figure, indicate that neither the Fisher 
et al. (2005) nor the Tout (1991) distribution properly matches 
the distribution we observe in GCs. 

5.2. The total binary fraction 

The procedure described in the previous section allowed us to 
directly measure the fraction of binaries with g>0.5 without 



any assumptions regarding/(q). On the other hand, because of 
the photometric errors, binaries with small mass ratios {q<0.5) 
are indistinguishable from single MS stars in this dataset, there- 
fore, any attempt to determine the total fraction of MS-MS bi- 
naries without assumption on the mass-ratio distribution is im- 
possible with this approach. 

The approach we follow to estimate the total fraction of bi- 
naries is similar to that used by Sollima et al. (2007) and con- 
sists of assuming a form for/(q). Since none of the two mass- 
ratio distributions available from literature properly matches 
the observed distribution in order to estimate the total frac- 
tion of binaries (/™ T ), we extrapolated the results of Sect. 15.11 
adopting a flat f(q) also for binary systems with ^<0.5; i. e. , 
we assumed a constant mass-ratio distribution for all q values. 
In this case as f^°' 5 = f^.' the total fraction of binaries is 
simply 

fTOT _ 2 o f^ >0 - 5 

•> bin ' ' .bin 

The final /™ T are listed in the fifth column of Table 2 for the 
rc, the rc-HM, the r oliM sample, and the WFC field. 
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Fig. 24. As in Fig.[23]for the remaining 30 GCs. 



For completeness, we note that, according to Fisher et 
al. (2005), 66.5% of binary systems have mass ratio larger 
than 0.5. Hence, assuming a Fisher et al. (2005) mass ration 
distribution, the total fraction of binaries should be: 

/i ot,f = i-504 /r 5 . 

If we assume that binary stars are formed by random associa- 
tions between stars of different masses, only 24% of binaries 
have <7>0.5, and the total fraction of binaries becomes: 
f^ = 4.167 



5.3. Comparison with previous measurements of the 
binary fraction in GCs 

To date, the fraction of binaries has been measured for 30 GCs. 
In Table 3 we list the photometric binary fraction in Galactic 
GCs from previous measurements and available in the litera- 
ture. Although for some GCs of our sample the fraction of bi- 
naries were already estimated in previous works, caution must 
be used to compare the results presented in this paper with lit- 
erature ones. In particular, it should be noted that the inferred 
values of the total fraction of binaries are tightly related to the 



assumed f(q). Many authors just determined lower limits to the 
binary fraction, as they studied binary systems with large q that 
are clearly separated from single MS stars. Without any indi- 
cation on the mass-ratio interval analyzed, a quantitative com- 
parison of results with these studies is not possible. 



From the comparison between Table 2 and 3 we note that 
in some cases the fraction of binaries measured in the same 
cluster region by different authors strongly differs from the re- 
sults presented here. As an example, in the case of NGC 6752, 
Rubenstein & Bailyn (1997) estimated an high fraction of bi- 
naries in the core (/t,i n =0.27+0.12), in disagreement with the 
results presented in this paper (/|™ T =0.03+0.03) and in Milone 

et al. (2010) 5 =0.03+0.01). To investigate these differ- 
ent results, Milone et al. (2010) re-examined the Rubenstein & 
Bailyn (1997) findings first by analyzing the same data with the 
improved photometric techniques that are now available, and 
then using the better datasets that have been collected more re- 
cently. They concluded that the disagreement comes from the 
use of the stellar photometry tools they used, which allow a 
better separation of stellar blends. Similarly, the large fraction 
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Table 2. Fraction of binaries with mass ratio q > 0.5, q > 0.6 and q > 0.7, and total fraction of binaries measured for rc, ?"c-hm, 
and r„HM sample and in the whole ACS/WFC field of view for the target GCs. We have analyzed only the ACS/WFC region with 
radial distance from the cluster center larger than a minimum radius /? mm whose value, in arcminutes, is listed in the first column. 



ID 


REGION 


^)>0.5 
■'bin 


^q>0.6 
-'bin 


M>01 
-*bin 


fI° T 

-'bin 


ARP2 

«mi„=0.0() 


r c sample 
rc-HM sample 
r oH M sample 
WFC field 


0.093+0.010 
0.119±0.023 

A AA 1 i A AO 1 

0.091 ±0.031 
0.096+0.009 


0.076±0.007 
0.093±0.017 
0.086±0.024 
0.079±0.006 


0.055±0.005 
0.056±0.012 
0.081±0.017 
0.057±0.004 


0.186±0.020 
0.238±0.046 
0.182±0.062 
0.192±0.018 


E3 

«min=0.00 


rc sample 
r C -HM sample 
r HM sample 
WFC field 


0.360±0.043 
0.317±0.203 
0.082±0.107 
0.347±0.041 


0.350±0.042 
0.147±0.171 
0.103±0.107 
0.336±0.039 


0.247±0.035 
0.264±0.171 
0.029±0.075 
0.237±0.033 


0.720±0.086 
0.634±0.406 
0.164±0.214 
0.694±0.082 


NGC 104 

« min =0.83 


r c sample 
rc-HM sample 
r oH M sample 

WFC field 


0.009±0.003 
0.009+0.003 


0.007±0.003 
0.007±0.003 


0.005±0.003 
0.005±0.003 


0.018±0.006 
0.018±0.006 


NGC 288 

«min=0.00 


rc sample 
r C -HM sample 
r oHM sample 
WFC field 


0.056+0.005 
0.054±0.007 

A AA^ i A AIM 

0.092+0.040 
0.056+0.004 


0.050±0.004 
0.045±0.005 

A AO 'I i A A 1 £1 

0.032±0.016 
0.048±0.003 


0.041±0.003 
0.030±0.004 

A AO 1 i A A 1 1 

0.02 1 ±0.01 1 
0.038±0.003 


0.112±0.010 
0.108±0.014 

A 1 O A i A AO A 

0.184±0.080 
0.112±0.008 


NGC 362 

«mi„=0.42 


r c sample 
rc-HM sample 
r oH M sample 

WFC field 


0.025±0.004 

A A 1 £1 i A AAO 

0.016±0.003 
0.020+0.003 


0.018±0.003 
0.01 1 ±0.003 
0.013±0.003 


0.010±0.003 

A AAO i A AA1 

0.008±0.003 
0.008±0.003 


0.050±0.008 

A AOO i A AA£" 

0.032±0.0()6 
0.040±0.006 


NGC 1261 

«mi„=0.00 


rc sample 
t-HM sample 
r oHM sample 

WFC field 


0.023±0.009 
0.032±0.004 

A AO. A i A A AO 

U.U20±U.UU3 
0.024±0.003 


0.023±0.006 
0.028±0.003 

A Al O i A /lA 1 ) 

0.016+0.003 
0.021 ±0.003 


0.021 ±0.005 
0.021±0.003 
0.012±0.003 
0.015±0.003 


0.046±0.018 
0.064±0.008 

A A>1A i A AA£" 

0.040 ±0.006 
0.048±0.006 


NGC 1851 
p — ft 


r c sample 
r C -HM sample 
r oH M sample 

wfL neia 


A AAO i A /lAO 

0.008+0.003 

U .UUo±U ,UU J 


A AAO i A AAO 

0.008 ±0.003 

U.UUo+U.UUj 


A AA£" i A AAO 

0.006 ±0.003 

U .UUO±U .UU D 


A A 1 £1 i A AA£" 

0.016±0.006 
u.uio±u.uuo 


NGC 2298 
« mi „=0.00 


rc sample 
rc-HM sample 
r oHM sample 
WFC field 


0.077±0.009 
0.056±0.007 
0.047±0.004 
0.073+0.004 


0.066±0.006 
0.047±0.005 
0.034±0.003 
0.054±0.003 


0.041±0.004 
0.036±0.004 
0.023±0.003 
0.036±0.003 


0.154±0.018 
0.112±0.014 
0.094±0.008 
0.146±0.008 


NGC 3201 

«min=0.00 


rc sample 
rc-HM sample 
r ohm sample 
WFC field 


0.064+0.004 
0.054±0.006 

0.061 ±0.003 


0.056±0.003 
0.039±0.004 

0.051 ±0.003 


0.042±0.003 
0.026±0.003 

0.037±0.003 


0.128±0.008 
0.108±0.012 

0.122±0.006 


NGC 4147 

«mi„=0.00 


r c sample 
rc-HM sample 
r HM sample 
WFC field 


0.131±0.047 
0.017±0.011 
0.019±0.006 
0.029±0.005 


0.103±0.036 
0.041±0.007 
0.019±0.003 
0.027±0.003 


0.044±0.021 
0.036±0.005 
0.012±0.003 
0.020±0.003 


0.262±0.094 
0.034±0.022 
0.038±0.012 
0.058±0.010 


NGC 4590 

«min=0.()0 


rc sample 
rc-HM sample 
r oH M sample 

WFC field 


0.057±0.006 
0.040±0.004 
0.053±0.007 
0.049±0.003 


0.054±0.004 
0.037±0.003 
0.038±0.005 
0.044±0.003 


0.040±0.003 
0.023±0.003 
0.025±0.003 
0.030±0.003 


0.114±0.012 
0.080±0.008 
0.106±0.014 
0.098±0.006 


NGC 4833 


r c sample 


0.033±0.004 


0.029±0.003 


0.021 ±0.003 


0.066±0.008 



r C -HM sample 
r oHM sample 
WFC field 



0.020±0.003 0.018±0.003 0.014±0.003 0.040±0.006 
0.029±0.003 0.025±0.003 0.018±0.003 0.058±0.006 



NGC 5024 



^min-0.75 



rc sample 
rc-HM sample 
r ohm sample 
WFC field 



0.028±0.003 
0.033±0.003 
0.031 ±0.003 



0.021 ±0.003 
0.024±0.003 
0.023±0.003 



0.014±0.003 
0.019±0.003 
0.017±0.003 



0.056±0.006 
0.066±0.006 
0.062±0.006 



NGC 5053 



« min =0.00 



rc sample 
r C -HM sample 
r oHM sample 
WFC field 



0.072±0.005 
0.093±0.020 



0.058±0.004 
0.072±0.013 



0.038±0.003 
0.050±0.010 



0.144±0.010 
0.186±0.040 



0.073±0.005 0.059±0.004 0.039±0.003 0.146±0.010 



NGC 5272 

«mi„=0.00 



r c sample 
rc-HM sample 
r HM sample 
WFC field 



0.027±0.007 
0.012±0.003 
0.019±0.003 
0.017±0.003 



0.031 ±0.004 
0.011±0.003 
0.015±0.003 
0.015±0.003 



0.024±0.003 
0.010±0.003 
0.012±0.003 
0.012±0.003 



0.054±0.014 
0.024±0.006 
0.038±0.006 
0.034±0.006 



NGC 5286 



^min-0.83 



rc sample 
r C -HM sample 
r oHM sample 
WFC field 



0.011 ±0.003 
0.009±0.003 



0.008±0.003 
0.006±0.003 



0.007±0.003 
0.006±0.003 



0.022±0.006 
0.018±0.006 



NGC 5466 



«„i„=0.00 



r c sample 
rc-HM sample 
r ohm sample 

WFC field 



0.071 ±0.004 
0.055±0.008 
0.016±0.035 
0.066±0.004 



0.058±0.003 
0.049±0.006 
0.022±0.024 
0.055±0.003 



0.041±0.003 
0.029±0.004 
0.009±0.010 
0.038±0.003 



0.142±0.008 
0.110±0.016 
0.032±0.070 
0.132±0.008 



NGC 5897 



«,„ in =0.00 



rc sample 
rc-HM sample 
r HM sample 
WFC field 



0.059±0.003 
0.025±0.017 



0.051 ±0.003 
0.012±0.011 



0.037±0.003 
0.008±0.008 



0.118±0.006 
0.050±0.034 



0.058±0.003 0.049±0.003 0.035±0.003 0.116±0.006 



NGC 5904 
fl mi „=0.67 



r c sample 
rc-HM sample 
r ohm sample 
WFC field 



0.012±0.003 
0.006±0.009 
0.011 ±0.003 



0.007±0.003 
0.003±0.004 
0.007±0.003 



0.005±0.003 
0.005±0.003 
0.005±0.003 



0.024±0.006 
0.012±0.018 
0.022±0.006 
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Table 2. Cont. 



ID 


REGION 


bin 


M>0.6 
-'bin 


,q>0.7 
Jbin 


j-TOT 
■Ann 


NGC 5927 
« min =0.00 


r c sample 

'C—HM otUiipiC 

r oHM sample 
WFC field 


0.052±0.009 

o rv26+n 003 

0.006±0.003 
0.016+0.003 


0.037+0.007 
016+0 003 
0.006±0.003 
0.012±0.003 


0.030±0.006 
014+0 003 
0.004±0.003 
0.010±0.003 


0.104±0.018 
052+0 006 
0.012±0.006 
0.032±0.006 


NGC 5986 
S mi „=0.83 


rc sample 

iTe-t t¥il a csi tnnlp 
' C-HM »tUiipiC 

r HM sample 
WFC field 


0.005+0.004 
0.006±0.003 
0.006±0.003 


0.003±0.003 
0.003±0.003 
0.003+0.003 


0.003±0.003 
0.001±0.003 
0.002±0.003 


0.010±0.008 
0.012±0.006 
0.012±0.006 


NGC 6093 


rc sample 

rr-~ Tin i;;imnlp 
'C—HM otuupic 

r oHM sample 
WFC field 


0.006±0.003 
0.006±0.003 


0.006±0.003 
0.006+0.003 


0.004±0.003 
0.004±0.003 


0.012±0.006 
0.012±0.006 


NGC 6101 

«mi„=0.00 


r c sample 

Yf-t t¥tl a mnl P 
' C-HM »tUiipiC 

r HM sample 
WFC field 


0.050+0.004 
0.042±0.004 
0.054±0.007 
0.048±0.003 


0.043±0.003 
0.040±0.003 
0.039±0.005 
0.041 ±0.003 


0.031 ±0.003 
0.026±0.003 
0.021 ±0.003 
0.028±0.003 


0.100±0.008 
0.084±0.008 
0.108±0.014 
0.096±0.006 


NGC 6121 

« min =0.00 


rc sample 
**c-hm sample 
r oHM sample 

WFC field 


0.074±0.007 
051+0 005 

0.061 ±0.004 


0.073±0.006 
047+0 004 

0.055±0.004 


0.052±0.005 
030+0 003 

0.039±0.003 


0.148±0.014 
107+0 010 

0.122±0.008 


NGC 6144 

« mi „=0.00 


r c sample 

' C-HM »tuiipic 

r oH M sample 
WFC field 


0.066±0.006 
039+0 005 
0.030±0.007 
0.048±0.003 


0.059±0.005 
079+0 004 
0.021±0.005 
0.040±0.003 


0.046±0.004 
017+0 003 
0.010±0.004 
0.028±0.003 


0.132±0.012 
078+0 010 
0.060±0.014 
0.096±0.006 


NGC 6171 

« min =0.00 


rc sample 

'C—HM o«UiipiC 

r ohm sample 
WFC field 


0.093+0.011 
0.046+0.003 

0.054±0.003 


0.071 ±0.008 
0.035±0.003 

0.042±0.003 


0.052±0.007 
0.027±0.003 

0.032±0.003 


0.186±0.022 
0.092±0.006 

0.108±0.006 


NGC 6205 

«mi„=0.00 


r c sample 
r C-HM sample 
r„HM sample 
WFC field 


0.005+0.003 
0.006±0.003 
0.012±0.003 
0.007+0.003 


0.010±0.003 
0.004±0.003 
0.006±0.003 
0.006±0.003 


0.007±0.003 
0.004±0.003 
0.004±0.003 
0.005±0.003 


0.010±0.006 
0.012±0.006 
0.024±0.006 
0.014±0.006 


NGC 6218 

«mi„=0.00 


r c sample 
rc-HM sample 
r oHM sample 
WFC field 


0.057+0.005 
0.032±0.003 
0.011±0.013 
0.037±0.003 


0.046±0.004 
0.025±0.003 
0.007±0.009 
0.030±0.003 


0.034±0.004 
0.019±0.003 
0.004±0.007 
0.023±0.003 


0.114±0.010 
0.064±0.006 
0.022±0.026 
0.074±0.006 


NGC 6254 
« min =0.00 


rc sample 
r C-HM sample 
r ohm sample 
WFC field 


0.039±0.004 
0.022±0.003 
0.027±0.007 
0.029±0.003 


0.032±0.003 
0.017±0.003 
0.018±0.005 
0.023±0.003 


0.023±0.003 
0.012±0.003 
0.012±0.003 
0.016±0.003 


0.078±0.008 
0.044±0.006 
0.054±0.014 
0.058±0.006 



NGC 6341 r c sample 



r c _HM sample 0.010±0.003 0.007±0.003 0.005±0.003 0.020±0.006 





r ohm sample 


0.009±0.003 


0.007±0.003 


0.004±0.003 


0.018±0.006 


«mi„=0.42 


WFC field 


0.011 ±0.003 


0.008±0.003 


0.006±0.003 


0.022±0.006 


NGC 6352 


rc sample 


0.092±0.008 


0.078±0.007 


0.054±0.005 


0.184±0.016 




t-HM sample 


0.053±0.005 


0.041±0.004 


0.034±0.003 


0.106±0.010 




r ohm sample 


0.039±0.017 


0.026±0.014 


0.015±0.011 


0.078±0.034 


«mi„=0.00 


WFC field 


0.069±0.004 


0.055±0.003 


0.042±0.003 


0.138±0.008 


NGC 6362 


r c sample 


0.060±0.004 


0.044±0.003 


0.034±0.003 


0.120±0.008 




r C -HM sample 


0.021 ±0.005 


0.020±0.004 


0.016±0.003 


0.042±0.010 




r oHM sample 


0.032±0.037 


0.023±0.026 


0.043±0.024 


0.064±0.074 


«mi„=0.00 


WFC field 


0.046±0.003 


0.037±0.003 


0.029±0.003 


0.092±0.006 


NGC 6366 


r c sample 


0.099±0.007 


0.082±0.006 


0.064±0.006 


0.198±0.014 




rc-HM sample 


0.057±0.015 


0.035±0.012 


0.042±0.012 


0.114±0.030 




r HM sample 










«mi„=0.00 


WFC field 


0.092±0.007 


0.074±0.006 


0.059±0.005 


0.184±0.014 


NGC 6388 


rc sample 












rc-HM sample 












r oHM sample 


0.004±0.004 


0.006±0.003 


0.003±0.003 


0.008±0.008 


^min = 0.83 


WFC field 


0.008±0.004 


0.006±0.003 


0.003±0.003 


0.016±0.008 


NGC 6397 


r c sample 


0.035±0.018 


0.037±0.015 


0.037±0.013 


0.070±0.036 




rc-HM sample 


0.012±0.003 


0.010±0.003 


0.005±0.003 


0.024±0.006 




r HM sample 


0.014±0.026 


0.005±0.003 


0.002±0.003 


0.028±0.052 


« mi „=0.00 


WFC field 


0.012±0.003 


0.0 11 ±0.003 


0.007±0.003 


0.024±0.006 


NGC 6441 


rc sample 












rc-HM sample 












r oHM sample 


0.010±0.005 


0.008±0.004 


0.006±0.003 


0.020±0.010 


«mi„=1.00 


WFC field 


0.010±0.005 


0.008±0.004 


0.006±0.003 


0.020±0.010 


NGC 6496 


r c sample 


0.089±0.006 


0.073±0.005 


0.051 ±0.004 


0.178±0.012 




r C-HM sample 


0.077±0.008 


0.053±0.007 


0.036±0.006 


0.154±0.016 




r oH M sample 


0.046±0.024 


0.021±0.018 


0.015±0.015 


0.092±0.048 


«mi„=0.00 


WFC field 


0.069±0.005 


0.049±0.004 


0.031 ±0.003 


0.138±0.010 


NGC 6535 


r c sample 


0.046±0.016 


0.027±0.012 


0.014±0.008 


0.092±0.032 




rc-HM sample 


0.026±0.013 


0.018±0.009 


0.018±0.009 


0.052±0.026 




r oHM sample 


0.028±0.010 


0.016±0.007 


0.012±0.006 


0.056±0.020 


«mi„=0.00 


WFC field 


0.033±0.009 


0.021 ±0.006 


0.014±0.005 


0.066±0.018 
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Table 2. Cont. 



ID 


REGION 


Al>0.5 
-•bin 


,q>0.6 
■'bin 


M>0J 
bin 


fTOT 
-'bin 


NGC 6541 
« min =0.42 


r c sample 
rc-HM sample 
r oHM sample 

WFC field 


0.014±0.003 
0.010±0.003 
0.010±0.003 


0.010±0.003 
O.O05±O.O03 
0.007±0.003 


0.005±0.003 
0.001 ±0.003 
0.003±0.003 


0.028±0.006 
0.020±0.006 
0.020±0.006 



NGC 6584 


rc sample 


0.045±0.006 


0.045±0.004 


0.034±0.003 


0.090±0.012 




r c-HM sample 


0.036±0.007 


0.025 ±0.005 


0.020±0.003 


0.072±0.014 




r HM sample 


0.025±0.003 


0.016±0.003 


O.O09±O.O03 


0.050±0.006 


«mi„=0.00 


WFC field 


0.039±0.003 


0.030±0.003 


0.021±0.003 


0.078±0.006 


NGC 6624 


rc sample 
















rc-HM sample 


0.013±0.004 


0.002±0.003 


0.001±0.003 


0.026±0.008 




r ohm sample 


0.013±0.005 


O.O18±O.O05 


0.010±0.004 


0.026±0.010 


«mi„=0.42 


WFC field 


0.011±0.004 


0.012±0.003 


0.007±0.003 


0.022±0.008 


NGC 6637 


r c sample 


0.062±0.010 


0.060±0.007 


0.057±0.006 


0.124±0.020 




rc-HM sample 


0.029±0.004 


0.028±0.003 


0.020±0.003 


0.058±0.008 




r HM sample 


0.0 13 ±0.003 


O.O08±O.O03 


0.005 ±0.003 


0.026±0.006 


«mi„=0.00 


WFC field 


0.030±0.003 


0.024±0.003 


0.019±0.003 


0.060±0.006 


NGC 6652 


rc sample 


0.172±0.055 


0.091±0.038 


0.059±0.029 


0.344±0.110 




rc-HM sample 


0.052±0.006 


0.032±0.004 


0.018±0.003 


0.104±0.012 




r oHM sample 


0.027±0.006 


0.021 ±0.005 


0.016±0.004 


0.054±0.012 


«mi„=0.00 


WFC field 


0.042±0.004 


0.027±0.003 


0.018±0.003 


0.084±0.008 


NGC 6656 


r c sample 


0.023±0.003 


0.018±0.003 


0.013±0.003 


0.046±0.006 




rc-HM sample 


0.020±0.003 


0.015±0.003 


0.010±0.003 


0.040±0.006 




r oHM sample 














« mi „=0.00 


WFC field 


0.022±0.003 


0.017±0.003 


0.012±0.003 


0.044±0.006 


NGC 6681 


rc sample 
















rc-HM sample 


0.026±0.005 


0.013±0.003 


0.006±0.003 


0.052±0.010 




r oH M sample 


0.005±0.004 


0.0 11 ±0.003 


0.010±0.003 


0.010±0.008 


fi min =0.10 


WFC field 


0.019±0.003 


0.013±0.003 


0.008±0.003 


0.038±0.006 


NGC 6723 


rc sample 


0.031 ±0.004 


0.025±0.003 


0.020±0.003 


0.062±0.008 




rc-HM sample 


0.013±0.003 


0.011±0.003 


0.006±0.003 


0.026±0.006 




r oH M sample 


0.017±0.004 


0.008±0.003 


0.008±0.003 


0.034±0.008 


«min=0.00 


WFC field 


0.023±0.003 


0.017±0.003 


0.013±0.003 


0.046±0.006 


NGC 6752 


r c sample 


0.017±0.016 


0.011 ±0.008 


0.006±0.004 


0.034±0.032 




rc-HM sample 


O.O05±O.O03 


0.004±0.003 


0.002±0.003 


0.009±0.006 




r HM sample 


— 


— 


— 


— 


*mm=0.00 


WFC field 


0.005±0.003 


0.004±0.003 


0.002±0.003 


0.010±0.006 


NGC 6779 


rc sample 


0.050±0.009 


0.050±0.006 


O.038±0.005 


0.100±0.018 




rc-HM sample 


0.028±0.003 


0.022±0.003 


0.017±0.003 


0.056±0.006 




r oHM sample 


0.023±0.003 


0.016±0.003 


0.012±0.003 


0.046±0.006 


«min=0.00 


WFC field 


0.028±0.003 


0.022±0.003 


0.017±0.003 


0.056±0.006 


NGC 6809 


r c sample 


0.040±0.003 


0.031 ±0.003 


0.023±0.003 


0.080±0.006 




rc-HM sample 


— 


— 


— 


— 




r oHM sample 


— 


— 


— 


— 


« min =0.00 


WFC field 


0.040±0.003 


0.031 ±0.003 


0.023±0.003 


0.080±0.006 


NGC 6838 


rc sample 


0.152±0.017 


0.120±0.015 


0.080±0.012 


0.304±0.034 




rc-HM sample 


0.110±0.008 


0.100±0.007 


0.072±0.006 


0.220±0.016 




r oHM sample 


0.104±0.014 


0.084±0.012 


0.076±0.011 


0.208±0.028 


«min=0.00 


WFC field 


0.117±0.007 


0.101±0.006 


0.074±0.005 


0.234±0.014 


NGC 6934 


r c sample 
















rc-HM sample 


0.032±0.003 


0.027±0.003 


O.O17±O.O03 


0.064±0.006 




r uHM sample 


0.020±0.004 


0.019±0.003 


0.012±0.003 


0.040±0.008 


«min=0.42 


WFC field 


0.023±0.003 


0.021±0.003 


0.013±0.003 


0.046±0.006 


NGC 6981 


r c sample 


0.049±0.009 


0.053±0.006 


0.041 ±0.005 


0.098±0.018 




rc-HM sample 


0.031 ±0.008 


0.035±0.006 


0.031 ±0.004 


0.062±0.016 




r HM sample 


0.034±0.006 


0.028±0.004 


0.019±0.003 


0.068±0.012 


«mi„=0.00 


WFC field 


0.038±0.004 


0.037±0.003 


0.029±0.003 


0.076±0.008 


NGC 7078 


rc sample 
















rc-HM sample 


0.010±0.005 


0.012±0.003 


0.009±0.003 


0.020±0.010 




r„ H M sample 


0.018±O.O03 


0.014±0.003 


0.007±0.003 


0.036±0.006 


^min=0-83 


WFC field 


0.017±0.003 


0.014±0.003 


0.007±0.003 


0.034±0.006 


NGC 7089 


r c sample 
















rc-HM sample 


0.032±0.006 


0.018±0.004 


O.O09±O.O03 


0.064±0.012 




r HM sample 


0.011±0.003 


0.009±0.003 


0.005±0.003 


0.022±0.006 


« mi „=0.83 


WFC field 


0.013±0.003 


0.009±0.003 


0.005±0.003 


0.026±0.006 


NGC 7099 


rc sample 


0.035±0.015 


0.033±0.015 


0.010±0.003 


0.070±0.030 




rc-HM sample 


0.012±0.003 


0.010±0.003 


0.008±0.003 


0.024±0.006 




r HM sample 


0.013±0.003 


0.009±0.003 


0.007±0.003 


0.026±0.006 


«min=0.00 


WFC field 


O.012±0.003 


0.010±0.003 


0.008±0.003 


0.024±0.006 


PALOMAR 1 


r c sample 


0.333±0.096 


0.311±0.092 


0.244±0.079 


0.666±0.192 




r c-HM sample 


0.130±0.042 


0.116±0.037 


0.093±0.033 


0.260±0.084 




r oH M sample 


0.095±0.031 


0.089±0.027 


0.070±0.023 


0.190±0.062 


« mi „=0.00 


WFC field 


0.146±0.027 


0.136±0.024 


0.108±0.021 


0.292±0.054 


PALOMAR 12 


r c sample 


0.130±0.057 


0.130±0.045 


0.104±0.037 


0.260±0.114 




rc-HM sample 


0.175±0.018 


0.144±0.015 


O.108±O.O13 


0.350±0.036 




r oHM sample 


0.066±0.019 


0.055±0.014 


0.044±0.012 


0.132±0.038 


«mi„=0.00 


WFC field 


0.137±0.013 


0.114±0.011 


0.087±0.009 


0.274±0.026 


TERZAN 7 


rc sample 


0.187±0.017 


0.159±0.013 


0.140±0.011 


0.374±0.034 




rc-HM sample 


0.084±0.016 


0.092±0.013 


0.073±0.010 


0.168±0.032 




r oH M sample 


0.088±0.011 


0.075±0.008 


0.051 ±0.006 


0.176±0.022 


« mi „=0.00 


WFC field 


0.117±0.008 


0.104±0.006 


O.083±0.005 


0.234±0.016 


TERZAN 8 


rc sample 


0.083±0.011 


0.072±0.008 


0.056±0.006 


0.166±0.022 




rc-HM sample 












r ohm sample 


0.059±0.009 


0.047±0.006 


O.O37±O.O05 


0.118±0.018 


«mi„=0.00 


WFC field 


0.067±0.007 


0.056±0.005 


0.044±0.004 


0.134±0.014 
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Fig. 25. Black points show the normalized fractions of bina- 
ries, Vbin measured in five mass-ratio intervals as a function of 
q for all the GCs studied in this paper. To compare the fraction 
of binaries in different clusters we divided Vb; n by two times the 
fraction of binaries with q>Q.5. For clarity, black points have 
been randomly scattered around the corresponding q value. Red 
points with error bars are the means in each mass-ratio bin, 
while the gray line is the best fitting line, whose slope is quoted 
in the inset. 




Fig. 27. Frequency of binaries divided by two times the fraction 
of binaries with q>Q.5 expected for the cases of Fisher et al. 
(2005) mass-ratio distribution (upper panel) and from the dis- 
tribution obtained by randomly extracting secondary stars from 
a Kroupa (2002) IMF (bottom panels). Red points with the er- 
ror bars are the mean values of the observed binary frequency 
normalized by two times / b q ; ^ and have been already plotted in 
Fig. [25] The slope of the best-fitting least-squares gray straight 
lines and the reduced-^ 2 obtained from the comparison of the 
observed and theoretical distribution are quoted in the figure. 



0.2 - 



~i — 1 — 1 — 1 — i — 1 — 1 — 1 — r 
Fisher et al. (2005) 



-| 1 1 1 1 1 1 1 1 1 1 r 

random extraction 

.from a Kroupa (2002) IMF 



0.47<M,<0.76 M„ 



0.2 



0.6 



0.8 



Fig. 26. Upper panel: Mass-ratio distribution derived by Fisher 
et al. (2005). Lower panel: Mass-ratio distribution simulated 
from random extraction from a Kroupa (2002) IMF in the pri- 
mary star mass intervals quoted in the inset. 



of binaries detected by Albrow et al. (2001) in NGC 104, and 
Fisher et al. (1995) are not confirmed by our study. 

Sollima et al. (2007) have recently measured the fraction of 
binaries in the core of 13 low-density GCs by using the same 



images studied in this paper. First, they analyzed the color dis- 
tribution of MS stars to directly derive the minimum fraction of 
binary systems required to reproduce the observed CMD mor- 
phologies, then they inferred two different estimates of the total 
fraction of binaries by assuming the mass-ratio distribution ob- 
tained from random extractions from a de Marchi et al. (2005) 
IMF, and from the distribution measured by Fisher et al. (2005). 

Even if we have shown that the Fisher et al. (2005) distribu- 
tion is not consistent with what found in the present work and 
because the images are the same as in this paper, for a meaning- 
ful comparison with Sollima et al. (2007), in Fig.|28]we com- 
pare the total fraction of binaries in the core that we obtained 
by assuming the Fisher et al. (2005) distribution (red circles) 
with the values from Sollima et al. (2007). Blue triangles cor- 
respond to the binary fraction estimated in this paper assuming 
a flat q distribution. For eight out of thirteen GCs, results are in 
agreement, at the level of less than three <x. In the cases of ARP 
2, NGC 6101, NGC 6723, NGC 6981, and Terzan 7 the frac- 
tion of binaries measured in this work is systematically smaller 
than those found by Sollima and collaborators. 

5.4. The binary fraction as a function of primary-star 
mass 

In this Section we investigate the distribution of binary sys- 
tems as a function of the magnitude. To do this, we calculated 
the fraction of binaries over the entire WFC/ACS field of view 
in the three magnitude intervals, containing all the single MS 
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Finally, we used isochrones to estimate the average mass of 
the single stars and the primary component of binary systems 
in the regions Ab, Aj, and Af. To do this we converted the mean 
F814W magnitudes of the single stars contained in each of 
these regions into masses through the Dotter et al. (2007) mass- 
luminosity relations. Figure[33l shows the ratio/^ ( ^ f //^ ' 5 as 
a function of the average mass estimated above and suggests 
that the binary fraction is nearly fiat in the analyzed mass inter- 
val. 

We recall here, that the results presented in this subsection 
come from the analysis of the binary fractions measured over 
the entire ACS/WFC field of view. Due to the relatively small 
numbers of binaries, we did not extended this analysis to each 
group of rc, the /"c-HMi and the r HM stars. 



Fig. 28. Comparison of the core binary fractions of 13 GCs 
measured in this paper (red circles) and in Sollima et al. (2007, 
black circles). In both cases has been assumed the mass-ratio 
distribution from Fisher et al. (2005). Blue triangles indicate 
the binary fraction estimated in this work assuming a flat q dis- 
tribution. 



stars and the binary systems with a primary star: [0.75,1.75], 
[1.75,2.75], [2.75,3.75], F814W magnitudes below the MSTO 
respectively. In the cases of NGC 6388 and NGC 6441 we 
used smaller magnitudes intervals of [0.75,1.25], [1.25,1.75], 
[1.75,2.25] F814W magnitudes below the MSTO. We divided 
the CMD regions A and B defined in Sect. [4] and illustrated in 
Fig. fTTI into three subregions (named Ab, A;, Af and B^, B u Bf) 
as shown in Fig. [29] and calculated the fraction of binaries in 
each magnitude interval see eq. Q] 

Results are shown in Figs. [30] and [3TI where we plot the 
fractions of binaries with mass ratio q>0.5 calculated in three 
magnitude bins as a function of the difference between the 
mean F814W magnitude of the bin and the F814W magnitude 
of the MS turn off (Ampsww)- Red points indicate the fraction 
of binaries in the full interval [0.75:3.75] ([0.75:2.15] for NGC 
6388 and NGC 6441), while the shadowed area indicates the 
error associated to this measure. 

In general we find no evidence for a significant trend in 
the fraction of binaries with magnitude (which is a proxy for 
primary mass), as suggested by the reduced-^ 2 values quoted in 
Figs. [30] and [3T| Montecarlo simulations show that in the case 
of a flat distribution the 50% and 99% of objects have^- 2 values 
smaller than 1.1 and 5.5 respectively. Possible exceptions to 
this rule of a flat trend are represented by NGC 5897 and NGC 
6652 for which we have estimated^- 2 values higher than 5.5. 
And large x 2 > 5.0 are obtained also for NGC 6144, NGC 
6637, and NGC 6723. 

In order to further analyze the general trend of the binary 
fraction with the magnitude for all the GCs studied in this paper 
we divided the values of/^ ' 5 measured in each magnitude bin 
by the fraction of binaries with q>Q.5 in the interval between 
0.75 and 3.75 magnitudes below the MS turn off. Results are 
shown in Fig. [32] and confirm that the fraction of binaries is 
nearly flat in the analyzed magnitude range. 




FB14W 

Fig. 32. Fraction of binaries with q>Q.5 measured in three 
magnitude intervals (black points) as a function of Ampsnw 
for the 59 GCs studied in this work. To compare the fraction 
of binaries in different clusters we have divided the fraction of 
binaries in each bin by the value of/2^ ' measured in the inter- 
val between 0.75 and 3.75 F814W magnitudes below the MS 
turn off. For clarity black points have been randomly scattered 
around the corresponding Ampsww value. Red points with er- 
ror bars are the average binary fractions in each interval while 
the gray line is the best fitting least-square line whose slope is 
quoted in the inset. 



5.5. The radial distribution 

In order to investigate how the fraction of high-mass-ratio bi- 
naries depends on the radial distance, we divided the ACS field 
of view into four concentric annuli, and calculated the fraction 
of binaries by following the recipes described in Sects. [4] and 
15.21 We chose the size of the annulus such that the number of 
stars that populate the CMD region A is equal in each of them. 
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Fig. 29. Dark and light gray areas indicate the CMD regions used to measure the fraction of binaries in three magnitude intervals. 
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Fig. 30. Fraction of binaries with mass ratio q>0.5 for 29 GCs measured in three magnitude intervals (black points) and in the 
interval between 0.75 and 3.75 F814W magnitudes below the MS turn off (red points) as a function of Arapsww- Horizontal 
segments indicate the magnitude coverage corresponding to each point. 



Results are shown in Figs. [34] and [35] where we plotted 
fS* ' 5 as a function of the explored radial distance for all the 
GCs studied in this paper, and confirm that, in most of the GCs 
where the fraction of binaries has been calculated both in the 



core and in the outer regions, binaries are significantly more 
centrally concentrated than single MS stars. 

In Fig. [36] for the 43 GCs studied in this paper for which 
we measured the core binary fraction, we plot the fraction of 
binaries with q>0.5 in units of core-binary fraction as a func- 
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Fig. 31. As in Fig.|30]for the remaining 30 GCs. In the cases of NGC 6388 and NGC 6441 the binary fraction has been measured 
in the interval between 0.75 and 2.25 F814W magnitudes below the MS turn off. 
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Fig. 33. Fraction of binaries with q>0.5 measured in three 
magnitude intervals and normalized by f£-' 5 (black points) as a 
function of the mass of the primary star for the 59 GCs studied 
in this work. 



Fig. 36. Binary fraction with q>0.5 (in units of core binaries) 
as a function of the distance from the cluster center in units of 
core radii 



tion of the radial distance in units of core radii While it would 
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Fig. 34. Fraction of binaries with q>0.5 as a function of the radial distance from the cluster center for 29 GCs. The dotted and 
dashed vertical lines mark the core and the half mass radius respectively. Black filled-circles show binary fractions in four radial 
intervals while red open-symbols indicate the binary fraction for the rc, rc-HM, and r HM sample. Horizontal black segments 
indicate the radial coverage corresponding to each point while observational errors are plotted as vertical lines and shadowed 
areas. 



be naive to assume a similar radial distribution of binaries for 
all Galactic GCs, it seems clear that the fraction of binaries 
normalized to the core binary fraction is correlated with radius, 
and that the binary fraction typically decreases by a factor of 
~2 at two core radii with respect to the core binary fraction. 
The latter behavior was also suggested by Sollima et al. (2007) 
on the basis of their analysis of the radial distribution of bina- 
ries in seven GCs. In the cluster envelope, the binary-fraction 
trend with radius tends to flatten. 

5.6. Correlation between the binary fraction and the 
parent-cluster parameters 

In this section we investigate whether the binary fraction is cor- 
related with any of the physical and morphological parameters 
of their host GCs. In particular, our analysis makes use of the 
following quantities: 

Relative ages. We used the most recent age measures by 
Marfn-Franch et al. (2009). Relative ages were obtained from 



the same photometric database used in this paper by compar- 
ing the relative position of the clusters' MS turnoffs, using MS 
fitting to cross-compare clusters within the sample. Typical er- 
rors on the relative age measurements are between 2 % and 7 
%. We also used absolute ages from Salaris & Weiss (2002) 
and De Angeli et al. (2005). Absolute ages are not available 
for 15 GCs, namely: E3, NGC 4147, NGC 4833, NGC 5024, 
NGC 5286, NGC 5927, NGC 5986, NGC 6144, NGC 6388, 
NGC 6441, NGC 6496, NGC 6541, NGC 7089, PAL 1, and 
TERZAN 8. 

Metallicity. We performed our analysis with both the metal- 
licity scales defined by Zinn & West (1984) and Carretta & 
Gratton (1997), which were also used by Marin-Franch et al. 
(2009) to determine relative ages. 

BSS Frequency. We used the counts of BSS derived by Moretti 
et al. (2008) from the WFPC2 photometric catalogs published 
by Piotto et al. (2002). In particular, we used the normalized 
number of BSS, which is the absolute number of BSS in a given 
region divided by the total luminosity coming from the stars in 
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Fig. 35. As in Fig.[34]for the remainin 



the same region (in unit of 10 4 L Q ). 
Rate of stellar collisions per year. 

King, Surdin & Rastorguev (2002) have shown that the 
rate of stellar collisions per cluster and per year is Tc = 
5 x lO^^rc) 1 / 2 , where So is the central surface brightness in 
units of L Q ypc~ 2 and rc is the core radius in units of parsecs. 
We calculated the probability (F*) that a given star will have 
a collision in 1 yr, by dividing the collision rate by the total 
number of stars in the cluster. This is calculated by assuming a 
mass-luminosity ratio of 2 and a mean mass for colliding stars 
of 0.4 M . 

We also compared the measured fraction of binaries with 
the encounter frequency adopted by Pooley & Hut (2006) in 
the approximation used for virialized systems: pQ 5 r c where rc 
is the core radius and po the central stellar density. 

The other parameters involved in this analysis are the ab- 
solute visual magnitude My, the ellipticity (e), the central con- 
centration (c), the core relaxation timescale, t c , the half-mass 
relaxation timescale Th m , and the logarithm of the central lumi- 
nosity density po, and are taken from the Harris (1996) com- 
pilation. We also used three different parameters related to the 



cluster HB morphologies, as discussed in Sect. 15.6.31 Ellipticity 
(e) measurements are not available for six clusters, namely 
ARP2, E3, NGC 288, PALOMAR 12, Terzan 7 and Terzan 8. 

Figures [3714471 show the monovariate correlations. Note 
that, in our study of the core population of binaries, we did not 
include the post-core-collapse (PCC) GCs, because, for these 
objects, the definition of core radius is not reliable (Trager et al. 
1993). Specifically, PCC clusters are marked with red crosses 
in these figures but are not used to study the statistical signif- 
icance of the correlations. Figures [3714471 show that there are 
no significant correlations between the binary fractions and the 
cluster ellipticity, core and half mass relaxation time, central 
concentration and metallicity as suggested by the small values 
of the Pearson correlation coefficient. Some marginal correla- 
tion with the central density can not be excluded. 

In the following we will discuss some of the relevant rela- 
tions between the cluster parameters listed above and the frac- 
tion of binaries calculated in three radial regions defined in 
Sect. [4] A noteworthy correlation of the binary fraction is with 
the central velocity dispersion (r — 0.6). as shown in Fig. [37ll38l 
and [39] The central velocity dispersions <x v come from Meylan 
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Fig. 37. Fraction of binaries with q > 0.5 in the core as a function of some parameters of their host GCs. Clockwise: ellipticity, 
central concentration, central velocity dispersion, logarithm of the central luminosity density, half-mass and core relaxation 
timescale, and metallicity. In each panel we quoted the Pearson correlation coefficient (r). PCC clusters are marked with red 
crosses and are not used to calculate r (see text for details). 
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Fig. 38. As in Fig. [37] for the rc-HM sample. 



(1989), and are available only for a subsample of the GCs that 
are studied in the present work. 

5. 6.1. /bin versus M v , T t , and BSS frequency 

The most significant correlation we found is the correlation be- 
tween the cluster binary fraction and its absolute magnitude. 
Clusters with fainter absolute luminosity (smaller mass) have 
higher binary fractions. This correlation is present in all binary 
groups, i. e. for binaries inside the cluster core, for binaries 
located between the core and the half-mass radius, for bina- 
ries outside the half-mass radius, and for binaries with q>0.5 
(Fig. l40b . An anticorrelation between the fraction of binaries 
and the mass of the host GCs is predicted by theoretical mod- 
els (Sollima 2008, see also Fregeau et al. 2009). These authors 



suggest that this correlation could be the due to the fact that 
cluster mass and the efficiency of binary destruction have the 
same dependence on the cluster density and velocity disper- 
sion. 

This anticorrelation might extend to open cluster masses. In 
fact, Sollima et al. (2010) found a dependence of the fraction 
of binaries and the cluster mass in a sample of five open clus- 
ters. Sollima et al. (2010) suggests that the binary disruption 
within the cluster core is the dominant process that determine 
the fraction of binaries in star clusters. 

Noteworthy, a similar anticorrelation between the fre- 
quency of BSSs and the absolute luminosity of the parent clus- 
ter has been found by Piotto et al. (2004), Leigh, Sills, & 
Knigge (2007), and Moretti et al. (2008). Interestingly enough, 
Fig- EH shows that the fraction of binaries is indeed correlated 
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Fig. 39. As in Fig. [37]for the r Q uM sample. 
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Fig. 40. Upper-left: Fraction of binaries with q > 0.5 in the core as a function of the absolute visual magnitude of the host GC. 
Dashed line is the best fitting straight line whose slope (s) and intercept (i) are quoted in the figure together with the Pearson 
correlation coefficient (r). PCC clusters are marked with red crosses and are not used to calculate neither the best-fitting line nor 
r. For completeness in the up per- right panels we show the same plot for the fraction of binaries with q > 0.6, and q > 0.7. Lower 
panels: Fraction of binaries with q > 0.5 in the rc-HM ( left) and r a uM (fight) sample as a function of My. 



with the fraction of BSSs. Sollima et al. (2008) observed a similar correlation between the BSS specific frequency and the 
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fraction of binaries in the core of 13 low-density Galactic GCs. 
These authors suggested that the evolution of primordial bina- 
ries could be the dominant BSS formation process (see also 
Knigge et al. 2009 and Leigh, Sills & Knigge 201 1). However, 
Davies et al. (2004) provided a simple model showing that the 
correlation between the BSS frequency and the cluster mass 
may be the result of the evolution of the binary fraction due 
to encounters. Here, we can only note that, figures. [42] and 
gammacolN seem to suggest a mild correlation between binary 
fraction and the collisional parameter, while there is no signifi- 
cant correlation between the BSS frequency and the collisional 
parameter (e. g. Piotto et al. 2004, Davies et al. 2004, Leigh 
Sills, & Knigge 2007, Moretti et al. 2008). It is clear that the 
connection between binaries and BSSs is far from trivial. The 
interpretation of the correlation of binary fraction with cluster 
parameters, and with BSS fraction is beyond the purposes of 
the present paper. 
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Fig. 41. Fraction of binaries with q > 0.5 as a function of 
the BSS frequency in the core. PCC GCs are marked with red 
points. 



sity clusters. In order to verify whether the binary fraction de- 
pendence on age is limited to low density clusters, in Fig. [46] 
we plot the binary fraction for the rc sample as a function 
of the age from De Angeli et al. (2005) and the relative age 
from Marin Franch et al. (2009) for clusters with central den- 
sity log(po) < 2.75 (same central density limit of Sollima et al. 
2007 sample). We also note that the youngest low density clus- 
ters in our sample have a larger binary fraction, but the fact that 
at least one old GCs (E3) hosts a large binary fraction suggests 
that more data are needed to confirm any systematic trend. 

5.6.3. HB morphology 

Binaries have been considered as a possible second parameter 
of the HB morphology by several authors. In particular, the ex- 
istence of a link between field B-type subdwarf (sdB) - which 
are the counterpart in the field of the extremely hot horizon- 
tal branch (EHB) stars in GCs - and binary systems is well- 
established, both on observational and theoretical grounds. A 
large population of binaries has been found among field sdBs 
(e. g. Napiwotzki et al. 2004 and references therein). However, 
the formation scenario of EHB stars in GCs may be different. 
In fact, several radial-velocity surveys for the measurement of 
the binary fraction among EHB stars have revealed a significant 
lack of binary systems (Moni Bidin et al. 2006, 2009). 

In order to investigate possible relations between the frac- 
tion of binaries and the HB shape we used three different pa- 
rameters: 

1) the median color difference between the HB and the RGB 
[A(V - /)], measured by Dotter et al. (2010) for 60 GCs using 
the same CMDs of this paper; 

2) The HB morphology index from Mackey & van den Bergh 
(2005); 

3) the effective temperature of the hottest HB stars (r e ff HB), 
measured by Recio-Blanco et al. (2006). A(V - /), HB index, 
and r e ff HB measures are available for 56, 55, and 28 GCs stud- 
ied in this paper. 

Monovariate relation between the fraction of binaries with 
q > 0.5 and these parameters are shown in Fig. [47] We find 
no significant relations between the fraction of binaries and the 
HB parameters, confirming the small or null impact of the bi- 
nary population on the HB morphology. Similar results are ob- 
tained for binaries with q > 0.6 and q > 0.7. 



5.6.2. / b i n versus age 

Figures [44] and [45] plot the fraction of binaries with q > 0.5 
in the core as a function of relative ages by Marin-Franch et al. 
(2009) and the absolute ages from and from Salaris & Weiss 
(2002) and De Angeli et al. (2005) respectively. There is no 
evident trend between ages and the binary fraction. 

Sollima et al. (2007) compared the fraction of binaries mea- 
sured in the core of thirteen clusters, with the cluster ages from 
Salaris & Weiss (2002) and De Angeli et al. (2005) and found 
an anticorrelation between age and binary fraction suggesting 
age as the dominant parameter that determines the fraction of 
binaries in a GC. Our sample of 59 GCs does not confirm such 
correlation. Sollima et al. (2007) sample is limited to low den- 



6. Summary 

In this paper we have analyzed the properties of the population 
of MS binaries of a sample of 59 GCs. The main dataset con- 
sists in the ACS/WFC images of the Globular Clusters Treasury 
project (GO10775, PI Sarajedini) that allowed us to obtain uni- 
form and deep photometry for an unprecedented number of 
GCs (see Sarajedini et al. 2007 and Anderson et al. 2008 for 
details). We have also used ACS/WFC, WFC3 and WFPC2 
data from the archive to obtain proper motions when images 
overlapping the GO10775 data are available. The CMDs have 
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Fig. 42. Fraction of binaries with q > 0.5,q > 0.6, and q > 0.7 in the rc region (upper panels) and fraction of binaries with 
q > 0.5 in the tc-hm and r D HM regions [bottom panels) as a function of the collisional parameter (F»). The adopted symbols are 
already defined in Fig.l40l 



been corrected for the effects of differential reddening and pho- 
tometric zero point variations due to small inaccuracies in the 
PSF model. 

We have measured the fraction of binaries with mass ratio 
<7>0.5 and estimated the total fraction of binaries for MS stars 
that are between 0.75 and 3.75 magnitudes fainter than the MS 
turn off. We have found that: 

- in nearly all the GCs the fraction of binaries is significantly 
smaller than in the field, where the binary fraction is larger 
than 0.5 (e. g. Duquennoy et al. 1991, Fisher & Marcy 
1992) with a few relevant exceptions (E3, Palomar 1) where 
the total binary fraction is greater than ~0.4. 

- We have obtained the fraction of binaries in five intervals 
of q (for q>0.5) and found that the mass-ratio distribution 
is generally flat. 

- There is no evidence for a significant correlation of the bi- 
nary fraction with primary mass of the binary system. 

- We measured the fraction of binaries in the cluster core, in 
the region between the core and the half-mass radius, and 
outside the half-mass radius and studied their radial distri- 
bution. Binary stars are more centrally concentrated than 
single MS stars with the fraction of binaries generally de- 



creasing by a factor of ~2 from the center to about two core 
radii. 

- We investigated monovariate relations between the fraction 
of binaries (in the rc, ?"c-hm, and r a uM sample) and the main 
parameters of their host GCs (absolute magnitude, HB mor- 
phology, age, ellipticity, metallicity, collisional parameter, 
half mass and core relaxation time, central concentration, 
central velocity dispersion, and central luminosity density). 

- We found a significant anticorrelation between the fraction 
of binaries in a GC and its absolute luminosity (mass). 

- We found a marginal correlation between the cluster central 
density and the central velocity dispersion. 

- We did not find any significant relation between the binary 
fraction and the HB morphological parameters. 

- We confirm a significant correlation between the fraction of 
binaries and the fraction of BSSs, indicating that the main 
formation mechanism of BSSs must be related to binary 
evolution. 
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Fig. 43. As in Fig. [42] In this case we used the encounter frequency adopted by Pooley & Hut (2006) in the approximation used 
for virialized systems. 
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APPENDIX A. Reliability of the measured binary frac- 
tion. 

In this appendix we investigate whether the fraction of bina- 
ries with q>0.5 that we measured with the procedure described 
in Sect. [4] are reliable or are affected by any systematic uncer- 
tainty due to the method we used. The basic idea of this test 
consists of simulating a number of CMDs with a given fraction 
of binaries, measuring the fraction of binaries in each of them, 
and comparing the added fraction of binaries with the measured 
ones. 

Simulation of the CMD. 

We started by using artificial stars to simulate a CMD made 
of single stars following the procedure already described in 
Sect. 14.21 To simulate binary stars to be added to the simulated 
CMD we adopted the following procedure: 

- We selected a fraction /j™ T of single stars equal to the frac- 
tion of binaries that we want to add to the CMD and de- 



rived their masses by using the Dotter et al. (2007) mass- 
luminosity relation. In our simulations we assumed the val- 
ues of/™ 1 =0.05, 0.10, 0.30, and 0.50. 

- For each of them, we calculated the mass M2 = q X Mi of 
the secondary star and obtained the corresponding rap8i4w 
magnitude. Its color was derived by the MSRL. For sim- 
plicity we assumed a flat mass-ratio distribution. 

- Finally, we summed up the F606W and F814W fluxes of 
the two components, calculated the corresponding magni- 
tudes, added the corresponding photometric error, and re- 
placed the original star in the CMD with this binary system. 



As an example, in the upper panels of Fig. [48] we show the 
artificial star CMD made of single stars only (left panel), and 
the CMD where we added a fraction /™ T =0.10 of binaries 
(right panel), for the case of NGC 2298. 
Simulation of the differential reddening. 
To probe how well the reddening correction works, we 
considered a simple model. The simulation of the differential 
reddening suffered by any single star is far from trivial as 
we have poor information on the structure of the interstellar 
medium between us and each GC. For simplicity, in this 
work we assumed that reddening variations are related to the 
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Fig. 44. Fraction of binaries with q> 0.5, q>0.6, and q>0.7 in the rc region (upper panels) and fraction of binaries with q>0.5 
in the rc-HM and r Q HM regions (bottom panels) as a function of the relative age measured by Marfn-Franch et al. (2009). The 
adopted symbols are already defined in Fig. [40] 



positions (X, Y) of each stars by the following relations: 

AE(B -V) = d(cos(X') + sin(Y')) 
where 

X = Citz(X - Xmax)/ (^max _ ^minX 
Y' = C2?r(Y - Fmax)/(^max - Imin)- 

Here Xmin.max and Imin.max are the minimum and the max- 
imum values of the coordinates X and Y, C\ is a free param- 
eter that determines the maximum amplitude E(B - V) varia- 
tion, and C2 governs the number of differential reddening peaks 
within the field of view. In this work, we used for each GC the 
value of C\ that ranges from 0.005 to 0.05 to account for the 
observed reddening variation in all the GCs, while we arbitrar- 
ily assumed three values of C2=3, 5, and 8 to reproduce three 
different fine-scales of differential reddening. As an example, 
in Fig. [49] we show the map of differential reddening added to 
the simulated CMD of NGC 2298 that is obtained by assuming 
Ci = 0.025 and C 2 = 5. 

Then, we have transformed the values of AE(B-V) cor- 
responding to the position of each stars into AAp606w, and 
AAp8i4w and added these absorption variations to the F606W 
and the F814W magnitudes. The CMD obtained after we 
added differential reddening is shown in the bottom left panel 



of Fig. [48] for NGC 2298. We applied to this simulated CMD 
the procedure to correct for differential reddening described in 
Sect. 13.11 and obtain the corrected CMD shown in the bottom 
right panel. For each of these binary-enhanced simulated 
CMD, we also generated a CMD made of artificial stars by 
following the approach described in Sect. 14.21 In our investi- 
gation we did not account for field stars. For each combination 
of the f£P T and C2 we have simulated 200 CMDs with random 
values of the C\. 

Measurements of the binary fraction. 

Finally, we used the procedure of Sect. [4] to measure the frac- 
tion of binaries with mass ratio g>0.5 defined as: 

/.q>0.5 _ ^SIMU _ ^ART 

■'bin N A N A 

"SIMU "ART 

where N^ lMV and N^ ml] are the numbers of stars in the re- 
gions A and B in the CMD, as defined in Fig.QT]in the binary- 
enhanced simulated CMD and A^ RT and A^ RT the numbers of 
stars in the same regions of the artificial stars CMD. 

Results are shown in Fig. l50l where we plotted the differ- 
ence between the measured and the input fraction of binaries 
versus the parameter C\ for four difference values of the input 
binary fraction. We found that, for input binary fraction of 0.05, 
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Fig. 45. As in Fig.[44]but in this case we used the age measures from Salaris & Weiss (2002) and De Angeli et al. (2005). 
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Fig. 46. Fraction of binaries with q> 0.5 in the rc sample for low density clusters (log(po) < 2.75) as a function of the relative 
age from Marin-Franch et al. (2009) (left panel) an absolute age from Salaris & Weiss (2002) and De Angeli et al. (2005) (right 
panel). 



0.10 and 0.30, the average difference are negligible (< 0.5%), and 8 ( indicated in Fig. [50] with red circles, gray triangles, 

as indicated by the the black lines and the numbers quoted in and black crosses, respectively) give indeed the same average 

the inset. In the case of a large binary fraction (/™ T =0.5) the differences. Our comparison between the fraction of binaries 

measured fraction of binaries with q>0.5 is systematically un- added to the simulated CMD and the measured ones demon- 

derestimated by ~0.03. Apparently our results do not depend strate that the fraction of binaries determined in this work and 
on the value of the parameter Ci_. Simulations with C2 =3, 5, 
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Fig. 47. Fraction of binaries with q> 0.5 as a function of the temperature of the hottest HB stars (bottom), the HB morphology 
index {middle), and the median color difference between the HB and the RGB {top). 



listed in Table 2 are not affected by any significant systematic 
errors related to the procedure we adopted. 

We have also determined the fraction of binaries in five 
mass-ratio intervals by following the approach described in 
Sect. 15.11 for real stars. To this aim, we have divided the re- 
gion B of the CMD defined in Sect. [TTJinto fi ve subregions as 
illustrated in Fig. [22] for real stars. The size of each region is 
chosen in such a way that each of them covers a portion of the 
CMD with almost the same area. The resulting mass-ratio dis- 
tribution is shown in Fig.[5T| where we have plotted the fraction 
of binaries per unit q as a function of the mass ratio. As already 
done in the case of real stars, to compare the mass-ratio dis- 
tribution in simulated CMDs with different binary fraction, we 
have divided Vbi„ by two times the measured fraction of bina- 
ries with q>Q.5. The best fitting gray line closely reproduce the 
fiat mass-ratio distribution in input with Vbi n =l- 

Finally we have measured in the simulated CMDs the frac- 
tion of binaries with q > 0.5 in three intervals [0.75,1.75], 
[1.75,2.75], and [2.75,3.75] F814W magnitudes below the 
MSTO. To do this we used the procedure already described 
in Sect. 15.41 for real stars and we have normalized the f^°' 5 
value measured in each magnitude bin by the fraction of bi- 
naries with q > 0.5 measured in the whole interval between 



0.75 and 3.75 F814W magnitudes below the MSTO. Results 
are shown in Fig. [52] where we have plotted the normalized bi- 
nary fractions as a function of Awipsnw- The best-fitting gray 
line is nearly flat, and well reproduces the input magnitude dis- 
tribution. 

These tests demonstrate that both the mass-ratio distribu- 
tion determined in Sect. l5.1l for the 59 GCs studied in this work 
and shown in Fig.[25]as well as the binary fractions measured in 
different magnitude intervals in Sect. 15. 41 are not biased by sig- 
nificant systematic errors related to the procedure we adopted. 
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ID 


/bin 


note 


region 


reference 


E3 


0.29±0.09 


lower limit 


within 2 core radius 


Veronesi et al. (1996) 


ARP2 


>0.08 


lower limit 


within core 


Sollima et al. (2007) 




0.329-0.521 


all 


within core 


Sollima et al. (2007) 


NGC 104 (47 Tucanae) 


0.14±0.04 


all 


within half-mass radius 


Albrow etal. (2001) 




>0.05 


lower limit 


outside half-mass radius 


De Marchi & Paresce (1995) 




-0.02 


all 


outside half-mass radius 


Anderson (1997) 


NGC 288 


>0.10 


lower limit 


1 -6 core radius 


Bolte (1992) 




0.10-0.20 


all 


within half-mass radius 


Bellazzini et al. (2002) 




A Al +0.1 
u ' ul -0.0 


all 


outside half-mass radiu 


Bellazzini et al. (2002) 




>0.06 


lower limit 


within core 


Sollima et al. (2007) 




0.116-0.145 


all 


within core 


Sollima et al. (2007) 


NGC 362 


0.21±0.06 


all 


within half-mass radius 


Fischer etal. (1993) 


NGC 2808 


0.20±0.04 




outside half-mass radius 


Alcainoetal. (1998) 




0.04±0.01 


all 


outside half-mass radius 


Milone et al. (2010) 


NGC 3201 


<0.10 


upper limit 


outside half-mass radius 


Cote et al. (2007) 


NGC 4590 


>0.09 


lower limit 


within core 


Sollima et al. (2007) 




0.142-0.186 


all 


within core 


Sollima et al. (2007) 


NGC 5053 


>0.08 


lower limit 


within core 


Sollima et al. (2007) 




0.110-0.125 


all 


within core 


Sollima et al. (2007) 


NGC 5139 (tj Centauri) 


<0.05 


upper limit 


outside half-mass radius 


Elson etal. (1995) 


NGC 5272 (M3) 


'low' 


all 


outside half-mass radius 


Gunn& Griffin (1979) 




-0.04 


all 


outside half-mass radius 


Pryor et al. (1988) 




0.14±0.08 


all 


inside half-mass radius 


Zhao & Bailyn (2005) 




0.02+°^ 


all 


outside half-mass radius 


Zhao & Bailyn (2005) 


NGC 5466 


>0.08 


lower limit 


within core 


Sollima et al. (2007) 




0.095-0.117 


all 


within core 


Sollima et al. (2007) 


NGC 5897 


>0.07 


lower limit 


within core 


Sollima et al. (2007) 




0.132-0.171 


all 


within core 


Sollima et al. (2007) 


NGC 6101 


>0.09 


lower limit 


within core 


Sollima et al. (2007) 




0.156-0.210 


all 


within core 


Sollima et al. (2007) 


NGC 6121 (M4) 


23+" J4 

u ' iJ -0.23 


all 


inside half-mass radius 


Cote & Fischer (1996) 




-0.02 


all 


outside half-mass radius 


Richer et al. (2004) 


NGC 6341 (M92) 


00 +U UJ 


lower limit 


outside half-mass radius 


Anderson (1997) 


NGC 6362 


>0.06 


lower limit 


within core 


Sollima et al. (2007) 




0.118-0.127 


all 


within core 


Sollima et al. (2007) 


NGC 6397 


<0.07 


upper limit 


within half-mass radius 


Cool & Bolton (2002) 




0.051 ±0.0 10 


all 


within half-mass radius 


Davis et al. (2008) 




0.012±0.004 


all 


1.3-2.8 half-mass radii 


Davis et al. (2008) 


NGC 6656 (M22) 


03 + " 16 


all 


outside half-mass radius 


Cote et al. (1996) 


NGC 6723 


>0.06 


lower limit 


within core 


Sollima et al. (2007) 




0.161-0.218 


all 


within core 


Sollima et al. (2007) 


NGC 6752 


0.27±0.12 


all 


within core 


Rubenstein & Bailyn (1997) 




0.03±0.01 


lower limit, q>0.5 


within core 


Milone et al. (2010) 




02 +0 - 16 

u ' -0.02 

0.01±0.01 


all 


between core and half-mass radius 


Rubenstein & Bailyn (1997) 




lower limit, q>0.5 


between core and half-mass radius 


Milone et al. (2010) 


NGC 6792 


'low' 


all 


outside half-mass radius 


Catelan et al. (2008) 


NGC 6809 (M55) 


>0.06 


lower limit 


within core 


Sollima et al. (2007) 




0.096-0.108 


all 


within core 


Sollima et al. (2007) 


NGC 6838 (M71) 




all 


within half-mass radius 


Yan & Mateo (1994) 


NGC 6981 


>0.10 


lower limit 


within core 


Sollima et al. (2007) 




0.281-0.399 


all 


within core 


Sollima et al. (2007) 


NGC 7078 (M15) 


-0.07 


all 


within half-mass radius 


Gebhardt etal. (1994) 


NGC 7099 (M30) 


<0.05 


upper limit 


outside half-mass radius 


Alcainoetal. (1998) 


PALOMAR 12 


>0.18 


lower limit 


within core 


Sollima et al. (2007) 




0.408-0.506 


all 


within core 


Sollima et al. (2007) 


PALOMAR 13 


>0.30±0.04 


lower limit 


inside -18 core radii 


Clark etal. (2004) 


TERZAN 7 


>0.21 


lower limit 


within core 


Sollima et al. (2007) 




0.509-0.649 


all 


within core 


Sollima et al. (2007) 



Table 3. Collection of literature binary fraction estimates. For each GC we listed the measured fraction of binaries (/bin), specified 
if the latter is a lower limit, an upper limit or a measure of the total fraction of binaries. We also indicate the spatial region where 
this measure was done and give the reference. 
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Fig. 48. Artificial stars CMD for NGC 2298 (upper-left) and 
simulated CMD with a fraction of /™ T 0.10 of binaries added 
(upper-right). Bottom panels show the simulated CMD after 
we added differential reddening (left) and the simulated CMD 
after the correction for differential reddening (right). 
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Fig. 50. Difference between the measured fraction of binaries 
and the fraction of binaries in input as a function of the param- 
eter C\ for four difference values of the input binary fraction. 
Black lines indicate the average difference. Red circles, gray 
triangles and black crosses indicate simulations with C2 =3, 5, 
and 8 respectively. 
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Fig. 49. Bottom-left: Map of differential reddening added to the 
simulated CMD of NGC 2298. The gray levels indicate the red- 
dening variations as indicated in the upper-right panel. Upper- 
left and bottom-right panels show AE(B - V) as a function of 
the Y and X coordinate respectively for stars into 8 vertical and 
horizontal slices. 
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Fig. 51. Fractions of binaries per unit q measured in five mass- 
ratio intervals as a function of q for all the simulated GCs. To 
compare the q distribution in simulated clusters with different 
fraction of binaries, we divided Vb; n by two times the fraction 
of binaries with q>Q.5. For clarity, black points have been ran- 
domly scattered around the corresponding q value. The means 
normalized binary fractions in each mass-ratio bin are repre- 
sented by red points with error bars, while the gray line is the 
best fitting line, whose slope is quoted in the inset. 
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Fig. 52. Fraction of binaries with <7>0.5 in three magnitude in- 
tervals as a function of Amp8i4w for all the simulated GCs. To 
compare the measured fraction of binaries in different clusters 
we have divided the measured binary fractions in each mag- 
nitude interval by the value of/j^ 5 measured in the interval 
between 0.75 and 3.75 F814W magnitudes below the MS turn 
off. Red points with error bars are the means normalized bi- 
nary fractions in each magnitude interval. The gray line with 
the quoted slope is the best-fitting least-squares line. 
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